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We study multiband semiconducting nanowires proximity-coupled with an s-wave superconductor 
and calculate the topological phase diagram as a function of the chemical potential and magnetic 
field. The non-trivial topological state corresponds to a superconducting phase supporting an odd 
number of pairs of Majorana modes localized at the ends of the wire, whereas the non-topological 
state corresponds to a superconducting phase with no Majoranas or with an even number of pairs of 
Majorana modes. Our key finding is that multiband occupancy not only lifts the stringent constraint 
of one- dimensionality, but also allows having higher carrier density in the nanowire. Consequently, 
multiband nanowires are better-suited for stabilizing the topological superconducting phase and for 
observing the Majorana physics. We present a detailed study of the parameter space for multiband 
semiconductor nanowires focusing on understanding the key experimental conditions required for the 
realization and detection of Majorana fermions in solid-state systems. We include various sources of 
disorder and characterize their effects on the stability of the topological phase. Finally, we calculate 
the local density of states as well as the differential tunneling conductance as functions of external 
parameters and predict the experimental signatures that would establish the existence of emergent 
Majorana zero-energy modes in solid-state systems. 



PACS numbers: 

I. INTRODUCTION 

The search for Majorana fermions has become an ac- 
tive and exciting pursuit in condensed matter physics [1 - 
[4]. Majorana fermions, particles which are their own 
antiparticles, were originally envisioned by E. Majorana 
in 1937 [5 in the context of particle physics (ie., the 
physics of neutrinos). However, the current search for 
Majorana particles is mostly taking place in condensed 
matter systems [6 , 7 where Majorana quasi-particles ap- 
pear in electronic systems as a result of fractionalization, 
and as emergent modes occupying non-local zero energy 
states. The non-locality of these modes provides the abil- 
ity to exchange and manipulate fractionalized quasiparti- 
cles and leads to non-Abelian braiding statistics [8,,14 . 
Hence, in addition to being of paramount importance 
for fundamental physics, this property of the Majoranas 
places them at the heart of topological quantum comput- 
ing schemes |T3llT5H29] . We mention that solid-state sys- 
tems, where the Majorana mode emerges as a zero energy 
state of an effective (but realistic) low-energy Hamilto- 
nian, enable the realization of the Majorana operator it- 
self, not just of the Majorana particle. Consequently, 
Majorana physics in solid-state systems is in fact much 
more subtle than originally envisioned by E. Majorana 
in 1937. For example, in condensed matter systems the 
non-local non-Abelian topological nature of the Majo- 
rana modes that are of interest to us is a purely emergent 
property. 

About ten years ago. Read and Green [9 discovered 
that Majorana zero-energy modes can appear quite nat- 
urally in 2D chiral p-wave superconductors where these 



quasiparticles, localized at the vortex cores, correspond 
to an equal superposition of a particle and a hole. A year 
later, Kitaev [11 introduced a very simple toy model for 
a ID Majorana quantum wire with localized Majorana 
zero-energy modes at the ends. Both these proposals 
involve spinless p-wave superconductors where one can 
explicitly demonstrate the existence of Majorana zero- 
energy modes by solving the corresponding mean field 
Hamiltonian. Recently, several groups [30l [3Tj suggested 
a way to engineer spinless p-wave superconductors in the 
laboratory using a combination of strong spin-orbit cou- 
pling and superconducting proximity effect, thus opening 
the possibility of realizing Majorana fermions in solid- 
state systems to the experimental field. The basic idea of 
the semiconductor /superconductor proposal [31 is that 
the interplay of spin-orbit interaction, s-wave supercon- 
ductivity and Zeeman spin splitting could, in principle, 
lead to a topological superconducting phase with local- 
ized zero-energy Majorana modes in the semiconductor. 
Since then, there have been many proposals for realizing 
solid-state Majoranas in various superconducting het- 
erostructures [3QH451. Among them, the most promis- 
ing ones involve quasi-lD semiconductor nanowires with 
strong spin-orbit interaction proximity-coupled with an 
s-wave superconductor [37, 38, 40 . The main advantage 
of this proposal is its simplicity: it does not require any 
specialized new materials but rather involves a conven- 
tional semiconductor with strong Rashba coupling such 
as In As or InSb, a conventional superconductor such as 
Al or Nb, and an in-plane magnetic field. High quality 
semiconductor nanowires can be epitaxially grown (see, 
for example, Ref. [46 for InAs and Ref. [47 for InSb) and 
are known to have a large spin-orbit interaction strength 
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a as well as large Lande ^f-factor (^inAs ~ 10 — 25 [48] 
and ^inSb ~ 20 — 70 [47 ). Furthermore, these materials 
are known to form interfaces that are highly transparent 
for electrons, allowing one to induce a large supercon- 
ducting gap A [49-51 . Thus, semiconductor nanowires 
show great promise for realizing and observing Major ana 
particles ^ iTj . It is important to emphasize that in the 
superconductor-semiconductor heterostructures the Ma- 
jorana mode is constructed or engineered to exist as a 
zero-energy state, and as such, it should be experimen- 
tally observable in the laboratory under the right condi- 
tions. 

In a strictly ID nanowire in contact with a super- 
conductor, the condition for driving the system into 
a topological superconducting phase [37l [38] is \Vx\ > 
Y^A^ + /i^, where Vx is the Zeeman splitting due to the 
in-plane magnetic field, /i is the chemical potential and A 
is the proximity-induced superconducting gap. Thus, the 
key to the experimental realization of Major ana fermions 
in this system is the ability to satisfy a certain set of re- 
quirements that ensures the stability of the Major ana 
bound states. The challenging task here is the ability 
to suppress effects of disorder, control chemical poten- 
tial fluctuations as well as fluctuations of other parame- 
ters. Given that realizing single channel (or one subband 
) nanowire is quite challenging, it is natural to consider 
semiconductor nanowires in the regime of multi-subband 
occupancy. It was shown in Ref. [40 that this is a promis- 
ing route and the existence of Majorana fermions does 
not require strict one dimensionality. In fact, the stabil- 
ity of the topological superconducting phase is enhanced 
in multiband nanowires due to the presence of "sweet 
spots" (multicritical points in the topological phase di- 
agram, see Ref. [52 for details) in the phase diagram 
where the system is most robust against chemical poten- 
tial fluctuations [40 . In this paper, we expand on these 
ideas and explore effects of various perturbations such as 
disorder in the superconductor and semiconductor, fluc- 
tuations in the tunneling matrix elements, etc. on the 
stability of the topological phase. Our goal is to iden- 
tify parameter regimes favorable for the exploration of 
the Majorana physics in the laboratory. We therefore 
use realistic physical models and realistic values of the 
parameters throughout this work so that our theoretical 
results are of direct relevance to experiments looking for 
Majorana modes in nanowires. 

The paper is organized as follows. We begin in Sec. [IT] 
by introducing a tight-binding model for the semicon- 
ductor nanowires, and derive the superconducting prox- 
imity effect. We show that electron tunneling between 
semiconductor and superconductor leads to important 
renormalization of the parameters in the semiconductor. 
In Sec. |III[ we calculate the low-energy spectrum in the 
regime of multiband occupancy and identify the topolog- 
ical phase diagram. In Sec. [IV| we study disorder effects 
on the stability of the topological phase. We consider 
several sources of disorder: short-range impurities in the 
superconductor, short-range and long-range impurities 



in the semiconductor as well as fluctuations in the tun- 
neling matrix elements across the interface. In Sec.[V] we 
present results for experimentally observable quantities 
(e.g., local density of states and tunneling conductance) 
calculated using realistic assumptions. Finally, we con- 



clude in Sec. VI with the summary of our main results. 



II. TIGHT-BINDING MODEL FOR 
SEMICONDUCTOR NANOWIRES 

A. Spin— orbit interaction and Zeeman terms 

Single-channel semiconductor (SM) nanowires have 
been recently proposed [37l [38] as a possible platform 
for realizing and observing Majorana physics in solid 
state systems. Obtaining strictly one-dimensional (ID) 
nanowires would raise significant practical challenges [49] 
^50|, but this requirement can be relaxed to a less strin- 
gent quasi-lD condition corresponding to multiband 
occupancy [40]. The physical system proposed for study- 
ing Majorana physics consists of a strongly spin-orbit in- 
teracting semiconductor, e.g., InAs and InSb, proximity- 
coupled to an s-wave superconductor (SC). The quasi-lD 
SM nanowire is strongly confined in the z direction, so 
that only the lowest corresponding sub-band is occupied, 
while the weaker confinement in the y direction is con- 
sistent with a few occupied sub-bands. Consequently, 
the linear dimensions of the rectangular nanowire satisfy 
the relation <^ Ly <^ which is the usual physical 
situation in realistic semiconductor nanowires. The low- 
energy physics of the SM nanowire is described by the 
Hamiltonian 

^nw = i^O + ^SOI = ^ tijcl^Cj^ - 4a^i^ 

Yl - + h.c] , (1) 
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where Hq includes the first two terms and describes hop- 
ping on a simple cubic lattice with lattice constant a and 
the last term represents the Rashba spin-orbit interaction 
(SOI). We include only nearest-neighbor hopping with 
tii-^s = —^0, where S are the nearest-neighbor position 
vectors. In Eq. (jlj), ct represents a spinor ct = (ct^,ct^) 

with ct^ being electron creation operators with spin a, 
fi is the chemical potential, a is the Rashba coupling 
constant, and & = (cr^^, a^y, a^) are Pauli matrices. In 
the long wavelength limit, ^ 0, our model reduces 
to an effective mass Hamiltonian with to = h^a~'^ /2rrf 
and Rashba spin-orbit coupling aR(kyCFx — ^x^y)^ where 
aR = aa. In the numerical calculations we use a set 
of parameters consistent with the properties of InAs, 
a = 5.3A, m* = 0.04mo, and aR = 0.1 eV-A. The posi- 
tion within the cubic lattice is described hy i = {ix, iy^iz) 
with 1 < ix{y,z) ^ ^x{y,z)- ^tic calculatious we have 
used = 10, Ny = 250, and A^^: between 10"^ and 2- 10^, 
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which corresponds to a nanowire with dimensions 1/^=5 
nm, Ly = 130 nm, and Lx between 5 /im and 10 /am. 

A brute force diagonahzation of Hamiltonian ([T]) on a 
lattice containing more than 10^ sites would be numeri- 
cally very expensive. More importantly, the relevant en- 
ergy scales in the problem are of the order of a few meV 
allowing to construct the low-energy model as we show 
below and significantly reduce the Hilbert space. The 
largest energy scale is given by the gap between the low- 
est sub-bands, which, for example, for the first and sec- 
ond subband and the parameters used in our calculation 
is AEsb ^ 1.6meV. Consequently, we are only interested 
in the low-energy eigenstates of the Hamiltonian ([T]). To 
obtain these states, we take advantage of the fact that 
the eigenproblem for Hq can be solved analytically and 
notice that i^soi can be treated as a small perturbation. 
Explicitly, the eigenstates of Hq are 



A=l 



sm 



Xa 



(2) 



where n = {nx^riy^nz) with 1 < nA < Nx^ and Xa is an 
eigenstate of the spin operator. The corresponding 
eigenvalues are 



-2to cos 



TTUx 

Nx^l 



AT 11 +CQS ^ ^ — 6 ] — /i, 

(3) 

where the chemical potential in the semiconductor ja is 
calculated from the bottom of the band. We project the 
quantum problem into the low-energy subspace spanned 
by the eigenstates ipria with energies below a certain cut- 
off value, Cn < Cmax, whcrc the cutoff energy emax is typ- 
ically of the order 15meV, i.e., one order of magnitude 
larger than the inter sub-band spacing. The number of 
states in this low-energy basis is of the order 10^ and thus 
the numerical complexity of the problem is significantly 
reduced. The matrix elements of the SOI Hamiltonian 



1 — 



Nx^l 



sm 



sm 



[x^y]}, (4) 



cos 



COS 



where the second term in the parentheses is obtained 
from the first term by exchanging the x and y indices. 

To realize nontrivial topological states that support 
Majorana modes, it is necessary that an odd number of 
sub-bands be occupied. [40] In the simplest case of a sin- 
gle sub-band nanowire, this can be achieved with the help 
of a Zeeman field {i.e. a spin splitting ) 



Zeeman 



(5) 



which opens a gap at small momenta and removes one of 
the helicities that characterize the spectrum of a Rashba- 
coupled electron system. In the multi-band nanowires 



the situation is more complicated (see below) but the 
Zeeman term is essential to avoid fermion doubling. The 
Zeeman term can be obtained using an external magnetic 
field applied along the x axis, T = gjiiBBx/2. When the 
chemical potential lies within one of the Zeeman gaps at 
= 0, the condition for odd sub-band occupation is sat- 
isfied and thus fermion doubling is avoided, which allows 
for the existence of Majorana modes at the ends of the 
nanowire. Note that in SM with a large g-factor, e.g., 
^inAs ^ 10 and ^inSb ~ 50, relatively small in-plane mag- 
netic fields can open a sizable gap without significantly 
perturbing superconductivity. This is a crucial ingredient 
of the present proposal which is particularly important 
in the context of the effect of disorder on the topologi- 
cal phase as we discuss in Sec. |IV[ For example, in In As 
a magnetic field Bx ~ IT corresponds to F ~ ImeV. 
Finally, we note that in the basis given by Eq. Q the 
Zeeman term has the simple form 

(Verier I ^Zeeman I V^riv) = T^w^aa', (6) 

where a = —a. 



B. Proximity— induced superconductivity 

In addition to spin-orbit interaction and Zeeman spin 
splitting, the only other physical ingredient necessary for 
creating the Majorana mode is the ordinary s-wave super- 
conductivity which can be induced in the semiconductor 
by the proximity effect, through coupling to an s-wave su- 
perconductor (SC). A model of the full system that sup- 
ports the Majorana modes contains, in addition to the 
nanowire and Zeeman terms, Eqns. ([T]) and ([5|, respec- 
tively, the Hamiltonian for the superconductor, Hsc^ and 
a term describing the nanowire-SC tunneling, i^nw-sc- 
We note that, to account quantitatively for the super- 
conductivity induced in the nanowire, one should also in- 
clude possible electron-phonon and electron-electron in- 
teractions within the SM itself, i^int- These interactions 
may enhance or inhibit the induced effect, depending on 
the details of the SM material [53l [54). In this paper 
we do not take into account effect of interactions in the 
semiconductor on the proximity effect and use a simple 
model for the proximity effect using tunneling Hamilto- 
nian approach [55 which is appropriate for the sample 
geometry considered here (thin semiconductor lying on 
top of the superconductor). The effects of interactions 
on the topological superconducting phase were recently 
considered in Refs. [52, 56-58 . In addition to affecting 
the proximity-induced SC gap, the repulsive Coulomb 
interactions among the SM electrons lead to an effective 
enhancement of the Zeeman splitting which might be fa- 
vorable for inducing topological superconductivity [58] , 

Thus, the total Hamiltonian for our model of semicon- 
ductor/superconductor heterostructure is given by 



■ HSC + ^nw-SC, (7) 
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where the tunnehng term reads 



-sc 



h.C 



(8) 



with Cicr and ajo- being electron destruction operators 
acting within the SM and^SC, respectively. We assume 
that the matrix elements Uj couple the sites of the SC 
located at the interface, j = (?"|| , ^^interface)/^^^ where ry 
is position vector in a plane parallel to the interface and 
^interface IS the Coordinate of the interface layer in a slab 
geometry, to the first layer of the semiconductor wire. 

The proximity effect can be now derived by integrating 
out superconducting degrees of freedom in Eq. ^ and 
considering the resulting effective low-energy theory for 
the SM. To identify the form of the effective low-energy 
Hamiltonian, we consider first the case of a single-band 
SM coupled to an s-wave SC through an infinite planar 
interface, then we address the specific issues related to 
multiband nanowires. 



2^1 - C/K with C = (A-£i7^)/A and being the Fermi 
level in the superconductor. With these approximations, 
the surface self-energy becomes 



E(cj) = \i\^v{eF) / deG{e,uj) 



(13) 
(14) 



In the homogeneous case the last term in Eq. (13) rep- 



resents a shift of the chemical potential and can be ne- 
glected as the chemical potential should be determined 
self-consistently by solving the appropriate equation for 
the fixed total electron density. 

We can now include the surface self-energy Ti{uj) into 
the SM Hamiltonian and study the effective low-energy 
model for the semiconductor. This can be done by inves- 
tigating the poles of the SM Green's function 



i. Infinite planar interface 

Within the tunneling Hamiltonian approach, the 
proximity-effect induced by an s-wave SC can be cap- 
tured by integrating out the superconducting degrees of 
freedom and calculating the surface self-energy due to 
the exchange of electrons between SC and SM [59j. We 
briefly review this approach here and use these results 
later when discussing the disordered s-wave superconduc- 
tors. The interface self-energy is given by 

^aa' {r, r', uj) = Tr^^,^ J(r, ri)G^^/ (n, ^2, a;)t'^(r2, r'), 

(9) 

where f(r, ri) is the matrix describing tunneling between 
semiconductor and superconductor. To illustrate the ba- 
sic physics, we use here the simplest form for tunnel- 
ing matrix elements t(r,ri) = i5{z)5{zi)S{r^^ — r[') with 
r'l and z denoting in-plane and out-of-plane coordinates. 
After some algebra, the surface self-energy Ti{r—r'^uj) is 
given by 



G-^{k,uj) =uj 1 



S(r-r^a;) = \i\^S{z)d{z' 



(27r)3 



)G(p,w), 
(10) 



and finally becomes in the momentum space 

np\^) = \i? j ^G{p,c.) (11) 

= \i\'f^de I '^5{e-i^)G{e,oj). (12) 

where A is half bandwidth. The density of states 
z/(5,pll) = / ^^{s — S,p) is usually a weakly- dependent 
function of momenta and energy u{eF,P^^) ^ ^{^f) = 



7A0 



(15) 
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where = h'^k'^ /2m^ — /i, aR is the Rashba coupling, 
r the strength of a Zeeman field oriented perpendicu- 
lar to the interface, Aq is the value of the supercon- 
ducting gap inside the SC, and 7 is the effective SM- 
SC coupling. In the calculations, in addition to the 
values specified in Sec. |II A[ i.e., m* = 0.04mo and 
aR = O.leV^A, we have Aq = ImeV. The effecting cou- 
pling 7 = t'^\tlj{iz = 1)\'^v{€f) depends on the trans- 
parency of the interface, t, the amplitude of the SM wave 
function at the interface, ip{iz = 1), and the local den- 
sity of states of the non-superconducting metal at the 
interface, which can be expressed in terms of the half- 
bandwidth A and the Fermi energy sf of the metal. [59] 
Note that the Green's function (15) is written in the 
Nambu spinor basis {u^^ u^^vi^ — '^^j^using the Pauli ma- 
trices Tx and ax that correspond to the Nambu and spin 
spaces, respectively. The identity matrices tq and ctq are 
omitted for simplicity. 

Explicit comparison between the effective theory de- 
scribed by Eq. (15) and microscopic tight-binding 



calculations [59 has shown remarkable agreement. A sim- 
ilar effective description has proven extremely accurate in 
describing the proximity effect induced at a topological 
insulator - superconductor interface. [60 Can the low- 
energy physics contained in Eq. (15) be captured by 
an effective Hamiltonian description? To address this 
question, we determine the poles of the Green's function 
at frequencies within the superconducting gap, < Aq, 
i.e., we solve the Bogoliubov-de Gennes (BdG) equation 
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det[G- 



-11 - 



0. Explicitly, we have 




Note that we have considered only the lowest energy 
mode. In Eq. (16) dynamical effects are generated 



by the frequency dependence of th e proximi ty-induced 
terms containing the expression 7/ y^Aq — cj^, with a rel- 
ative magnitude that depends on the SM-SC coupling 
strength. In general, we distinguish a weak-coupling 
regime characterized by 7 <: Aq and a strong coupling 
regime, 7 ^ Aq. In the weak coupling regime we expect 
negligible dynamical effects at all energies that are not 
very close to the gap edge, u = Aq. Neglecting the fre- 
quency dependence in the proximity-induced terms, the 
solution of Eq. ( 15 ) becomes 



(17) 

where Z=(l+7/Ao) ^<lis the quiasiparticle residue 
at zero energy. 

To evaluate the dynamical effects, we compare the 
BdG spectrum given by Eq. (17) with the full solution of 
Eq. (16). The results are shown in Fig. [l]for a weak cou- 
pling regime characterized by 7 = 0.25meV (top panel) 
and at large coupling, 7 = 2meV (bottom panel). Notice 
that, even for effective couplings larger than the gap, e.g., 
7 = 2Ao, Eq. (17) represents a very good approximation 



of the low energy spectrum. What about couplings that 
are much larger than Aq? To answer this question, let 
us consider the dependence of the low-energy spectrum 
on the coupling constant 7, and the Zeeman field P. For 
F = the spectrum is gaped and the minimum of the 
gap is located at /c = (see Fig. [T]). Applying a Zeeman 
field reduces this minimum gap continuously and at the 
critical value Fc = \/7^ + /i^ the spectrum becomes gap- 
less (see Fig. [T]). For F > Fc a gap opens again with a 
minimum at /c = for F > Fc and at a finite wave vector 
for large values of the Zeeman field. The dependence of 
the minimum gap on the Zeeman field is shown in Fig. 
[2] for three different values of the SM-SC coupling. The 
vanishing of the gap at the critical point F = Fc marks a 
quantum phase transition from a normal SC at low Zee- 
man fields to a topological SC ( when F > Fc).[371 [59] 
The change in the location of the quasiparticle gap from 
A: = to a finite wave vector is marked in Fig. |2] by a 
discontinuity in the slope at F* > Fc. The optimal value 
of the excitation gap in the topological phase is obtained 
for F F*. This optimal value depends weakly on the ef- 
fective coupling 7, but varies strongly with the spin-orbit 
coupling. 

From the above analysis we conclude that the strong- 
coupling regime characterized by 7 ^ Aq is not ex- 
perimentally desirable, as it would require extremely 
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FIG. 1: (Color online) Top: Low-energy BdG spectrum of 
a semiconductor with proximity-induced superconductivity. 
The effective SM-SC coupling is 7 = 0.25meV and the chem- 
ical potential is /i = 0. The induced ga p vanish es at /c = 
in the presence of a Zeeman field F = v^7^ + /i^ = 0.25meV. 
The full lines are obtained by solving Eq. |T6| , while the 
symbols represent the spectrum given by Eq. (|17|). Bottom: 
Low-energy BdG spectrum for 7 = 2meV, /x = and three 
different values of the Zeeman field. The filled area corre- 
sponding to energies above Ao = ImeV represents th e SC 
continuum. Note that for energies Ek < Ao/2, Eq. (17) rep- 
resents a very good approximation of the BdG spectrum even 
in the strong-coupling regime. 



high magnetic fields to reach the topologically nontrivial 
phase, i.e., F > \J ^ 11^ ^ Aq. In addition, it would 
be difficult to tune the chemical potential and drive the 
system into a topological superconducting phase for a 
large SM-SC coupling. Also, as follows from Eq. (17), the 



quasiparticle excitation spectrum decays with increasing 
7/A0 which leads to the reduced stability of the topo- 
logical phase against thermal fluctuations. Hence, an ex- 
perimentally useful interface should be characterized by 
an effective coupling 7 of order Aq or less, i.e., in the 
intermediate to weak-coupling regime. As shown above, 
in these regimes the BdG spectrum is accurately approx- 
imated by Eq. (17). Consequently, we can model the 



low-energy spin-orbit coupled semiconductor with prox- 
imity induced superconductivity using an effective tight- 
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FIG. 2: (Color online) Dependence of the minimum quasi- 
particle excitation gap in the BdG spectrum given by Eq. 
|T6| on the Zeeman field F for different SM-SC couplings. 
The chemical potential is /x = and the Rashba coefficient is 
= 0.15eV-Afor the curve represented by small (green) cir- 
cles and ar = O.leV-Afor the other th ree curves. The system 
becomes gapless at Fc = y^7^~+7?. The superconducting 
state with F < Fc is topologically trivial, while for F > Fc 
one has a topological superconductor that supports Majorana 
bound states. Note that the optimal quasiparticle gap for the 
topological SC has a weak dependence on the SM-SC cou- 
pling but varies strongly with the strength of the spin-orbit 
coupling. 



binding model given by the Hamiltonian 



M 



H' 



Zeeman 



Ha, 



(18) 



where the semiconductor term Hsu has the same form as 
the Hamiltonian i^nw foi" the nanowire given by Eq. ([T]) 
but with Nx oo and Ny oo and the Zeeman term 
^Zeeman IS givcu by Eq. (pi with ax cf^. In addition, 
all the energy scales are renormalized by a factor Z = (1+ 
7/Ao)~"^ < 1, i.e., to Zto, a Za, etc. The physical 
meaning of the factor Z written as Z = 7~^/(7~^ + 
A^^) is intuitively clear - it corresponds to a probability 
to find an electron in the semiconductor. The induced 
superconductivity is induced by the effective pairing term 



E (a4 



h.c. 



(19) 



with an effective SC order parameter A = 7Ao/(7 + Ao). 
With these choices, i^eff given by Eq. (18) has the same 



low energy spectrum as the one described by Eq. (17) 



2. The multiband case 

To account for the specific aspects that characterize 
the proximity effect in finite-size systems, we return to 
the details of deriving the effective low-energy Green's 



function description, e.g., Eq. (15), starting with the 
microscopic Hamiltonian After integrating out the 
SC degrees of freedom, the effective SM Green's function 
acquires a self-energy 

l)t(i,,i,) (20) 

*cc I'f'y ijx i3y 

X GsC , ^ 1 1 , ^interface )t{jx , jy ) '4^n' {jx JyA), 

where i/jniix^'^yA) are the orbital components of the 
eigenstates described by Eq. ([2| and Gsc{^, '^\\ i ^interface) 



3,0 is ^he SC Green's function, both evaluated at the inter- 



face. In Eq. (20) r\\ = {i 



- jy) and t(ix,iy) are 

matrix elements that couple two sites with in-plane co- 
ordinates {ix, iy) across the interface. In the discussion of 
the planar interface we implicitly assumed translational 
invariance for the SM — SC coupling, i.e., t{ix,iy) = t. 
Here, we consider position-dependent couplings and ar- 
gue that engineering interfaces with a transparency that 
varies across the wire, i.e., t{ix, iy) = t{iy)^ generates off- 
diagonal components of the effective SC order parameter 
that help stabilize the Majorana modes. [40] We note that 
variations of the coupling matrix element along the wire, 
i.e., in the x direction, act as an effective disorder poten- 
tial. We will address this issue below in Sec lIV CI 

The SC Green's function integrated over momenta can 
be written explicitly as (see Eq. (13)) 



SC 



-u{eF) 



C 



i-C^ 



(21) 



where ( = {A — £:f)/A with A being the half-bandwidth, 
Ef being the Fermi energy and Aq = ImeV being the 
s-wave SC gap. In the numerical calculations we have 
Ef — A/2, i.e., C = 0.5. Using the expression of the SC 
Green's function given by Eq. (21), the self-energy (20) 
becomes 



Ssc(n,nO = -7n^n; 



(22) 

with the implicit assumption that the wire is very thin. 
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1. Tfie coupling matrix in Eq. 

(23) 



1 



sm 



1 



{N. + 1)A 



(24) 



We note that for position-independent SM-SC couplings, 
t{iy) = t, the matrix 7 is proportional to the unit ma- 
trix and the effective low-energy Hamiltonian can be ob- 
tained along the lines of Sec. |IIB 1| However, for non- 
uniform couplings, 7n n' acquires off-diagonal elements 
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FIG. 3: (Color online) Upper panel: Non-uniform SM-SC 
coupling t{y)_ = ^P^y)] with a smooth profile p{y) and 

e = [t(0) - t(L^)]/t(0) = 0.8. Lower panel: Dependence of 
the induced gap on the non-uniformity of the coupling across 
the wire. The inhomogeneous proximity effect induces inter- 
sub-band pairing with An^n^ii ~ for > 0.7. 



that generates normal and anomalous inter-sub-band 
terms in the effective Hamiltonian via the self-energy 
(22). The relative magnitude of the off-diagonal terms 



depends on the non-homogeneity of the SM-SC coupling. 
To quantify this property, we consider a fixed profile 
p{y) with the property p{0) = and p{Ly) = 1 and 
the position-dependent tunneling t{y) = to[l — Op{y)]^ 
where < ^ < 1 is a parameter that measures the de- 
gree of non-uniforniity of the coupling. Shown in Fig 
(upper panel) is t{y) for = 0.8. As we show be- 
low, the non-uniform coupling induces an effective pair- 
ing An^n; = {ny\j Ao/{j + Ao)|np, where 7 is given by 
Eq. (24). The dependence of ^nyn'y on the coupling 
asymmetry parameter is shown in Fig [s] (lower panel). 
For uniform tunneling {0 = 0) the off-diagonal elements 
vanish and the diagonal elements become equal. By con- 
trast, strongly non-homogeneous tunneling (^ ^ 1) gen- 
erate off-diagonal terms A^^n' that reach about 50% 
of the diagonal contributions for neighboring sub-bands. 



n'y = n^y ± 1, and are much smaller for |nj^ — n^yl > 1. 

The low-energy effective Hamiltonian for the nanowire 
can b e derived following the scheme described in Sec. 
II Bl As before, at low energies {uo <C Aq) we can 



neglect the frequency dep endence of the dynamically- 
generated terms, i.e., 1/y^Aq — cj^ ^ I/Aq. However, 
due to the inter-sub-band coupling induced by non- 
homogeneous proximity effect, one cannot simply renor- 
malize the energy by a factor Z . This is due to the 
fact that the Green's function for the proximity-coupled 

nanowire, {G~^)rtn' = ^ ^nn' - (^nw + ^Zeeman)w " 

T^Q^ciji^n')^ contains a frequency-dependent term uj(l + 
^riyn'y/ ^q) that is not proportional to the unit matrix. 
However, we notice that we can define a matrix Z^/^ 
with the property Z^I^Q-^Z^I^ = uj - H^^. Explicitly, 
we have 




(25) 



Because det[Z^/2] > 0, 



{Z^^'^)nyn'y - 

where ^{iy) is given 
v/2/(7V^ + lsin[i^n^7r/(7V^ 
the renormalized Green's function satisfies the same BdG 
equation as the original Green's function, i.e., det[co' — 
^eff] = det[G~^] = 0. We conclude that the low-energy 
physics of a nanowire proximity-coupled to an s-wave SC 
can be described by an effective Hamiltonian i^eff that 
can be conveniently characterized by its matrix elements 
in the Nambu basis \\n(j)) = {ipncr^ —(^^^lia)^ provided 
by the eigenstates of Hq given by Eq. ([2|. Considering 
only the lowest band, i.e., n = {fix^riy^ 1), we can write 
explicitly 



{{na\\H,f,\\n'a'))= (26) 



where a 



(n^,r7i 



,,1), = (n^,m2y/, 1), and 
summation over the repeating indices m^, rriy is implied. 
In Eq. (26) the matrix Z^/^ is given by Eq. (25), the 
Hamiltonian for the nanowire is i^nw = + ^sob and 
the matrix elements of i^o, ^soi, and i^zeeman are given 
by equations ([3|, (|4]), and ([6|, respectively. We note that 
for a homogeneous SM-SC interface the normal contri- 
bution proportional to ^nyn'y becomes diagonal and can 
be absorbed in the chemical potential, but in general 
it generates inter-sub-band mixing. These induced off- 
diagonal terms can be significant in the strong-coupling 
limit. Finally, the effective SC order parameter is 



7A0 



Ao + 7 



(27) 



where 7 is given by Eq. 
significance as in Eq. (25) 



(24) and \ny) has the same 



Eq. (26 



The effective low-energy BdG Hamiltonian given by 
is the main result of this section. In the remain- 
der of this work we will study the low-energy physics of 
a nanowire with proximity-induced superconductivity by 
diagonalizing numerically this effective Hamiltonian ( 26 ) . 
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III. LOW-ENERGY SPECTRUM, MAJORANA 
BOUND STATES, AND PHASE DIAGRAM 

A. General properties of the BdG spectrum 



The effective BdG Hamiltonian ( 26 ) can be written as 



eff 



H' 



Zeeman 



(28) 



where i^nw and i^zeeman are renormahzed nanowire and 
Zeeman Hamiltonians, respectively, and H/^ is the ef- 
fective pairing with matrix elements —dn^n'^^aa' ^Uyu' • 
We note that the normal contribution proportional to 
^Uyfi'y from Eq. (26) is included in Hnw As mentioned 
above, for homogeneous SM-SC coupling the only effect 
of this term is to generate an overall shift of the en- 
ergy. For convenience, we eliminate this shift by adding 
a term Sn n'^^crcr'^C/\/l ~ to the effective Hamil- 

tonian (26), where 



E 



Ao 



(29) 



is the "average" effective pairing. Note that ^UyUy ^ 
for Uy ^ 1 and the diagonal contribution to the energy- 
shifting term is partially canceled even in the non- 
homogeneous case. Finally, to be able to compare re- 
sults corresponding to various degrees of inhomogeneity 
in the coupling, i.e., different 6 parameters, we define the 
average coupling strength as 



(30) 



To obtain better insight into the properties of the BdG 
Hamiltonian, we start with a non-superconducting sys- 
tem described by the Hamiltonian H^w + ^Zeeman and 
Lx ^ oo, i.e., an infinitely long renormalized nanowire 
placed into an effective magnetic field. The spectrum of 
the renormalized wire is shown in Fig. [4] for a chemi- 
cal potential /J./Ea = 5 and a Zeeman field T/Ea ~ 15, 
where Ea = m*aR ^ 0.6K is the characteristic spin- 
orbit coupling energy. Here and below we systematically 
use Ea as energy unit. The bare nanowire spectrum is 
renormalized due to a weak inhomogeneous coupling with 
a profile given in Fig. [S) l9 = 0.8, and 7 = 0.25Ao. For 
F = the sub-bands with a given value of Uy are dou- 
ble degenerate at /c^^ = 0, but this degeneracy is removed 
by the Zeeman field. However, for special values of F, 
sub-bands corresponding to different values of Uy may 
become degenerate at kx = (see Fig. |4|. If the chemi- 
cal potential has a value such that the degeneracy point 
occurs at zero energy, (F, /i) represent a so-called "sweet 
spot" . [40 Adding superconductivity, will now open a gap 
in the spectrum near E = 0. In the weak-coupling limit, 
one can determine the topological nature of the induced 




-0.04 -0.02 



FIG. 4: (Color online) Spectrum of an infinite non- 
superconducting wire in the presence of a Zeeman field. 
The energies are obtained by diagonalizing the Hamiltonian 
ifnw + i^Zeeman for /i — 5Ea and F ^ 15Ea and are mea- 
sured relative to the chemical potential. Due to the presence 
of the Zeeman field, each Uy band is split into two sub-bands 
marked {riy—} and {ny-\-}. Note that, for this value of F, the 
sub-bands {!+} and {2—} are degenerate at /c^c = 0. 



superconductivity by simply counting the number of sub- 
bands crossed by the chemical potential: an odd number 
corresponds to a topologically non-trivial SC, while an 
even number results in a standard superconductor. [40] 
Within this simplified picture, the "sweet spots" repre- 
sent critical points. As will be shown below, the prop- 
erties of the system in the interesting parameter regimes 
near the "sweet spots" are determined by the effective 
inter-band pairing, i.e., by the non-homogeneity of the 
SM-SC proximity effect. 

Next, we diagonalize numerically the full effective 
Hamiltonian (26) for a finite wire using the same set of 

5 and T/E^ 



control parameters, /J./Ea = 5 and T/E^ ~ 15. The cor- 
responding low-energy spectrum is shown in the upper 
panel of Fig. [5] The eigenstates are labeled by an inte- 
ger number n that has the same sign as the corresponding 
eigenvalue En. The spectrum is characterized by a gap 
A* ^ l.SEa and a pair of zero-energy Majorana bound 
states. To prove the localized nature of these states, we 
calculate the wave function amplitude and show that the 
Majorana modes are localized near the ends of the wire 
(Fig. [5] lower panel). By contrast, finite energy states ex- 
tend over the entire system. The oscillations of the wave 
function amplitudes is associated with the Fermi momen- 
tum kp ^ 0.02/a, as shown in Fig. [4] We note that 
amplitudes shown in Fig. [5] represent the particle com- 
ponent of the BdG wave- functions, i.e., |i/ntP + I'^nj,^- 
For a finite-energy state, e.g., n = 2 the corresponding 
total spectral weight is about 1/2, with the other half 
coming from from n = —2. The Majorana modes have 
a total weight of one, which corresponds to one physi- 
cal particle, but this weight is spatially separated into 
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two contributions localized near the ends of the wire. 
Removing the Majorana pair would require overlapping 
the two components, which cannot be done by local per- 
turbations. This is, of course, the topological immunity 
of the Majorana modes, which is crucial for topological 
quantum computation. The characteristic length scale 
for the localized modes is controlled by the minimum 
value of the quasi-particle gap in a wire with no ends, 
e.g., with periodic boundary conditions, A^. As will be 



shown below (see Sec IIIC), in a finite wire it is possible 
that bound states localized near the ends of the system 
have energies within the gap. When in-gap states are 
present, the lowest-energy localized state sets the value 
of the mini-gap. A* < A^. We emphasize that for a 
set of control parameters (F, /i) corresponding to a non- 
vanishing minimum quasi-particle gap A^, the mini-gap 
A* is always nonzero. The characteristic length scale for 
the localized modes diverges in the limit A* 0. Hence, 
the topological phase is protected as long as the quasi- 
particle gap remains finite. Consequently, to determine 
the stability of the Majorana bound states, our key task 
is to determine the dependence of A^ on various physical 
parameters and experimentally-relevant perturbations, 
e.g., chemical potential, Zeeeman field, SM-SC coupling, 
and charged impurity and coupling-induced disorder. 

At this point in our analysis it is important to clar- 
ify the role played by the parameters that incorporate 
the SM-SC proximity effect into the low-energy effec- 
tive theory. In particular, we address the following ques- 
tion: how does the low-energy spectrum depend on the 
strength of the SM-SC coupling (i.e., on 7), on the non- 
homogeneity of the coupling (6>), and on the dynamical 
effects included in the effective description (Z^/^)? The 
non-homogeneity of the coupling is responsible for gen- 
erating inter sub-band pairing in Eq. ( [27| ). These off- 
diagonal contributions play a minor role away from the 
"sweet spots", where they generate a small quantitative 
change of the quasi-particle gap. However, in the vicinity 
of the "sweet spots" A^ vanishes in the absence of inter 
sub-band pairing and the non-homogeneity of the cou- 
pling (i.e., ^ > 0) becomes crucial. As expected, inclu- 
sion of dynamical effects through Z^/^ renormalizes the 
energy scales. Without these effects, the mini-gap would 
increase monotonically with the coupling strength, but 
dynamical effects limit its maximum value. The optimal 
quasi-particle gap obtains in the intermediate coupling 
regime 7 ^ Aq . Further increase of the coupling leads to 
a decrease of the gap. To illustrate the features described 
above, we show in Fig. [6] low-energy spectra in the vicin- 
ity of the "sweet spot" (/i = 14.5^;^ F = 15.3Ea) at 
weak coupling (7 = 0.25Ao, top panel) and intermediate 
coupling (7 = Ao, bottom panel). Note that in the ab- 
sence of inter-band pairing, i.e., for homogeneous SM- 
SC coupling, the gap near the "sweet spot" collapses. 
Also, inclusion of dynamical effects at intermediate and 
strong coupling is the key for obtaining the correct en- 
ergy scales. Finally, the location of the "sweet spots" in 
the F — /i plane depends on the coupling strength and. 
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FIG. 5: (Color online) Upper panel: BdG energy spectrum 
for a finite superconducting wire obtained by numerical di- 
agonalization of i^eff- The parameters fi and F are the same 
as in Fig. |4]and n labels the eigenvalues of Hes staring with 
the lowest energy and has the same sign as En. The in-gap 
states are Majorana zero-energy modes. Lower panel: The 
particle-component of the wave-function amplitude for the 
lowest energy states. The Majorana modes (n = 1) are lo- 
calized at the ends of the wire, while the finite energy states 
extend over the entire wire. 



more generally, the location of phase boundaries depends 
on the strength of the SC proximity effect. 



B. Phase diagram for mult i— band superconducting 
nanowires 

Topological superconductivity and, implicitly, the Ma- 
jorana bound states are protected by the quasi-particle 
gap A^, as discussed above. The vanishing of sig- 
nals a transition between topologically nontrivial and 
topologically trivial superconductivity. In a multi-band 
system, such transitions can be caused, for example, by 
varying the Zeeman field while maintaining a fixed value 
of the chemical potential. The vanishing of A^ at cer- 
tain specific values of F reveals a sequence of alternating 
SC phases with trivial and nontrivial topological proper- 
ties. A natural question is whether different topologically 
nontrivial (or trivial) phases have exactly the same low- 
energy properties. While topologically identical, these 
phases may have have some distinct features, at least in 
certain parameter regimes. 

To address this question, we calculate the low-energy 
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FIG. 6: (Color online) Upper panel: Low-energy spectra for 
weak SM-SC coupling (7 = 0.25Ao) in the vicinity of the 
"sweet spot" (/x = 14. 5^^, F = 15.3Ea). The yellow circles 
correspond to inhomogeneous coupling with ^ = 0.8 and takes 
into account dynamical effects, the black squares are obtained 
by neglecting dynamical effects, i.e., (Z^^^)^^^/^ = ^^y^y' 
the red diamonds are for a homogeneous coupling. Lower 
panel: Same as in the upper panel for 7 = Ao near the "sweet 
spot" (/i = 14.5^c., F = 24.2^a). Note that inclusion of dy- 
namical effects through Z^^^ renormalizes the energy scales, 
while inhomogeneous coupling play a critical role in establish- 
ing a finite gap near the "sweet spot" . The location of a given 
"sweet spot" depends on the coupling strength. 



spectrum for a system with a fixed chemical potential, 
/i = 30£^Q,, and different values of F, one from each of 
four intervals separated by points characterized by a van- 
ishing gap. The results are shown in Fig. jT] The dis- 
tinctive feature of the spectra in Fig. [7| is their num- 
ber of zero-energy modes. In the presence of a very 
weak Zeeman field, the SC phase has trivial topology and 
no zero-energy modes. Increasing F generates a transi- 
tion to a topological SC phase characterized by one pair 
of Majorana bound states. The further increase of the 
applied Zeeman field produces alternating topologically 
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FIG. 7: (Color online) Sequence of low-energy spectra ob- 
tained for four different values of the Zeeman field separated 
by points with a vanishing mini-gap. The spectra with an odd 
number of pairs of zero-energy modes (N) characterize topo- 
logical SC phases, while those with an even N correspond to 
trivial SC phases. Note the overall decrease of the mini-gap 
with the Zeeman field. The system is characterized by the 
following parameters: /x = SOEa, 7 = 0.25Ao, and = 0.8. 



trivial and nontrivial SC phases with increasing number 
of zero-energy bound states. Topological SC phases are 
characterized by an odd number of pairs of Majorana 
bound states, while trivial SC phases have an even num- 
ber of pairs. 

This type of sequence of alternating SC phases is in- 
dependent of the chemical potential or the strength of 
the SM-SC coupling. As discussed above, in the weak- 
coupling limit 7^0 this behavior can be directly related 
to the number of sub-bands of an infinite wire that cross 
the chemical potential: an odd (even) number corre- 
sponds to a nontrivial (trivial) topological superconduc- 
tor. At the phase boundary between two SC phases with 
different topologies, two sub-bands become degenerate 
dit kx = 0. In addition, at certain special values of the 
control parameters F and /i two different phase bound- 
aries intersect leading to multicritical points. We will call 
these special crossing points x-points. The "sweet spots" 
mentioned above are examples of such x-points. The key 
question is how the phase boundaries evolve when we 
turn on the SM-SC coupling and, in particular, what the 
physics is in the vicinity of the x-points, see also Ref. [52] , 
To obtain the global phase diagram in the T — jn plane. 




FIG. 8: (Color online) Phase diagram of the multiband 
nanowire as function of the Zeeman field F and the chemical 
potential /i. The quasi-particle gap of the effective low- 
energy Hamiltonian ( |26| ) vanishes at the phase boundaries. 
Superconducting phases characterized by an odd (even) num- 
ber of pairs of zero-energy Majorana modes are topologically 
nontrivial (trivial). The coupling of the SM nanowire to the 
s-wave SC is characterized by 7 = 0.25Ao and = 0.8. The 
inhomogeneous coupling induces off-diagonal pairing A^^^/^ 
which removes the x-points creating regions with stable non- 
trivial (near the "sweet spots") or trivial SC. 



FIG. 9: (Color online) Phase diagram of the multiband 
nanowire at intermediate coupling. The SM-SC coupling is 
7 = Ao, while the other parameters are the same as in Fig. 
[S] Note that the "sweet spots" of the topological phase char- 
acterized by = 1 are significantly expanded as compared 
with the weak-coupling case shown in Fig. [8] As a result, 
tuning the Zeeman field in the vicinity of F = 30^a allows 
for huge variations of the chemical potential without crossing 
a phase boundary, i.e., without closing the gap. 



we determine the parameters that satisfy the condition 
A^(r,/i) = 0, i.e., we identify the phase boundaries for 
transitions between topologically trivial and non-trivial 
phases by imposing the condition of a vanishing quasi- 
particle gap. We show below that this approach is consis- 
tent with the calculation of the topological index M. (Ma- 
jorana number) which distinguishes topologically trivial 
and non-trivial SC phases [TT1|40]. The results are shown 
in Figjsjfor weak coupling (7 = 0.25Ao) and in Figj9]for 
intermediate coupling (7 = Aq). 

For 7 = O.25A0 (see Fig. [8|, the only significant differ- 
ence as compared to the weak-coupling picture presented 
above is the disappearance of the phase boundary cross- 
ings at the x-points. Instead, the region in the vicinity 
of these points is occupied by the phase that is robust 
against variations of the chemical potential. Near the 
"sweet spots" , this phase is the nontrivial topological SC. 
We note that the disappearance of phase boundary cross- 
ings is a direct result of the off-diagonal pairing induced 
by a non-uniform SM-SC coupling. Uniform tunneling 
{0 = 0) does not eliminate the x-points, independent of 
the coupling strength, but pushes them to higher values 
of the Zeeman field as 7 increases. Also, we note that the 
characteristic width of the stable phase in a given avoided 



crossing region is controlled by a specific off-diagonal 
component Any^n'^- For example, the "sweet spots" in- 
side the phase characterized by = 1 (see Fig. [8| are 
controlled by the dominant matrix elements A^^^^+i (see 
Fig. |3|, while the avoided crossings within the N = 2 
topologically trivial phase are controlled by matrix ele- 
ments Anyny-\-2^ that are typically smaller. Hence, we 
expect the strongest effect within the topological phase 
characterized by one pair of Majorana bound states. As 
this phase also requires relatively low Zeeman fields, it 
is the experimentally relevant phase for realizing and ob- 
serving Majorana fermions. At intermediate couplings, 
the N = 1 phase is pushed to slightly higher fields (see 
Fig. [9|. However, as a result of the effective phase space 
of the "sweet spots" expanding significantly, this regime 
presents the remarkable possibility of being able to vary 
the chemical potential over energy scales of the order 
lOmeV without crossing a phase boundary. As we will 
show below, this feature has major experimental implica- 
tions in the sense that the elusive Majorana mode is most 
likely to be experimentally realized in the laboratory in 
this particular physically realistic parameter regime. We 
note that this interesting parameter regime exists for the 
multiband situation. 
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FIG. 10: (Color online) Dependence of the gap on the Zee- 
man field for a fixed chemical potential. The upper panel 
corresponds to weak SM-SC coupling with 7 = 0.25Ao and 
6 — 0.8, while the lower panel is obtained for an intermediate 
coupling with 7 = Ao and = 0.8. Within the red (darker 
gray) regions superconductivity is topologically trivial, while 
in the yellow (light gray) regions the system is topologically 
nontrivial. These curves correspond to horizontal cuts in the 
phase diagrams shown in figures [S] and [9] respectively, through 
a sweet spot of the phase N = 1. Note that in the limit ^ ^ 0, 
i.e., for uniform SM-SC coupling, the width of the lower field 
topologically nontrivial region shrinks to zero for both cou- 
pling strengths. 



C. Dependence of the gap on the Zeeman field and 
the chemical potential 

The phase diagram provides information about the 
topological nature of the phase characterized by given 
sets of control parameters (F,/i). However, we are ulti- 
mately interested in the stability of the Majorana bound 
states that occur at the ends of a nanowire within the 
topologically nontrivial phase. As the Majorana zero- 
energy modes are protected by the quasi-particle gap, 
knowing the size of the gap and the dependence of 
on the control parameters is critical. To address this is- 
sue, we determine the dependence of the gap on both 
the Zeeman field at fixed chemical potential and on fi at 

In 



fixed F. The results are shown in figures [TO] and 11 



general, the gap is non-zero everywhere except at points 
corresponding to phase boundaries. In the vicinity of a 
point (Fc,/ic) with A*(Fc,/ic) = the minimum gap in 
an infinitely long wire occurs at /c^^ = 0. The dependence 
of this minimum on the Zeeman field is approximately 
linear in the deviation from the critical field, |F — F^ 



FIG. 11: (Color online) Dependence of the quasi-particle gap 
on the chemical potential for a fixed Zeeman field. The up- 
per panel corresponds to weak coupling, 7 = 0.25Ao, and 
F = 25Ea, while the lower panel is obtained at intermediate 
coupling, 7 = Ao, and F = 32Ea. The curves represent ver- 
tical cuts through the N = 1 topological phase in the phase 
diagrams shown in figures [S] and [5] At intermediate coupling 
(lower panel) the gap is finite over the entire chemical po- 
tential range, making the topological SC phase in a system 
with average chemical potential /2 ~ 60^^ is robust against 
variations of the chemical potential of the order 6/j = 5meV. 



shown in Fig. [2] Note that outside this linear regime, 
the minimum gap occurs at finite wave vectors. Finally, 
we note that, at intermediate couplings, the Zeeman field 
can be tuned so that the gap remains finite over a large 
range of chemical potentials (Fig. 



11 , bottom panel). 



(see Fig. 10). This generalizes the single-band results 



Such regimes are extremely stable against fluctuations of 
the chemical potential produced by disorder or other per- 
turbations, as we will show explicitly in the next section. 
We emphasize that the critical ingredients for realizing 
this regime are: i) the off-diagonal pairing obtained by a 
non-uniform SM-SC coupling, and ii) an effective aver- 
age coupling 7 of the order of the bare SC gap Aq. Note 
that the coupling strength 7 ex P/AL^ can be controlled 
by either modifying the tunneling or by changing the 
width of the nanowire in the direction perpendicular 
to the interface with the superconductor. 

As we mentioned above, the value of the gap in a fi- 
nite system is, in general, smaller that the minimum gap 
in a system with the same parameters but no ends, e.g., 
with periodic boundary conditions. We emphasize that 
this is not a finite size effect, but is due to the appear- 
ance of in-gap states that are localized near the ends 
of the wire. In the topologically non-trivial phase, the 
characteristic length scale for these states is the same 
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FIG. 12: (Color online) Top panel: Spectrum of a system 
characterized by /i = 54.5£^q-, F — 35^^, 7 — Aq, and ^ — 0.8. 
The states with n — ±1 are Majorana zero-modes, n — 2,3 
correspond to states localized near the ends of the wire, and 
states with n > 4 are extended states. The mini-gap is A* = 
E2 and the minimum quasi-particle (QP) gap is A^ = E4 
The finite-energy localized states can be viewed as precursors 
of the extra-pair of Majorana bound states that obtains for 
F > . Bottom panel: Dependence of the quasi-particle 

(QP) gap A^ and mini-gap A* on the Zeeman field for a 
system with the same parameters as in Fig. |1Q[ Note that 
A* = A^ in the vicinity of phase transition points. The 
quasi-particle gap A^ has the same values as in Fig. [To] 



as that of the Majorana zero modes and, for wires with 
lengths larger than this scale, their energy is indepen- 
dent of Lj;. To illustrate this behavior, we show in Fig. 
12 (top panel) the spectrum of a system characterized by 
= 54.bEa, r 35^;^,, 7 Ao, and 6 0.8. We no- 
tice a pair of zero-energy Majorana modes and a number 
of states with almost the same energy that are extended 
over the entire length of the wire, as we checked explic- 
itly. In addition, we notice pair of states with energy 
within the bulk gap. An analysis of the position depen- 
dence of |?/^(x)p reveals that these states are localized 
near the ends of the wire. The minimum energy of the 
extended states is equal to the quasi-particle gap A^, 
while the lowest energy of the bound states is equal to 
the mini-gap A*. We can understand these bound states 
as precursors of the extra-pair of Majorana zero-modes 
that characterize the topologically trivial phase that ob- 
tains for F > 4:0Ea- Increasing the magnetic field pushes 
down the energy of the bound states until it vanishes at 
the transition, when A* = A^ = 0. On the other side 



of the transition, an extra-pair of localized states will be 
characterized by zero energy (see Fig. [7|. The evolu- 
tion of the mini-gap with the Zeeman field is shown in 
the lower panel of Fig. 12 Note that in the vicinity of 
the transition points A* = A^, while deep inside the 
topological phase A* < A^. Similar behavior can be 
observed throughout the phase diagram, including the 
topologically trivial phases. 



D. Phase diagram for an effective three— band 
model using the topological invariant. 

In this section we consider an effective three-band toy 
model which allows one to qualitatively understand sev- 
eral features observed in the detailed numerical simula- 
tions discussed in the previous sections. 

The topological phase diagram for the multiband 
nanowire can be obtained analytically using topological 
arguments due to Kitaev [11 . Namely, the supercon- 
ducting phase hosting Majorana fermions has an odd 
fermion parity whereas the non-topological phase has 
even fermion parity. Thus, the phase diagram can be 
obtained by calculating the Z2 topological index M (Ma- 
jorana number) defined as 

M = sgn[Pf5(0)]sgn[Pf5(7r/a)] = ±1. (31) 

The change of A4 signals the transition between trivial 
{A4 = 1) and non-trivial phases {A4 = —1). The anti- 
symmetric matrix B in Eq.(31) represents the Hamilto- 



nian of the system in the Majorana basis [n and can be 
constructed by the by the virtue of particle-hole symme- 
try [37, 61 . Specifically, the particle-hole symmetry of 
the BdG Hamiltonian is defined as 



(32) 



where O = /7i^ is an anti-unitary operator with U and 
K representing unitary transformation and complex con- 
jugation, respectively. One can check that the matrix 
B{P) = HBdG{P)U calculated at the particle-hole invari- 
ant points where i^BdG(^) = ^BdG(— ^) is indeed anti- 
symmetric {P) = —B{P). In ID there are two such 
points: P = 0,^ with ^ being the momentum at the 
end of the Brillouin zone and a being the lattice spacing. 
The function Pf in Eq. ( [31] ) denotes Pfaffian of the anti- 
symmetric matrix B. In the continuum approximation, 
where the lattice spacing a ^ and P = 7r/a^oo, the 
value of sgn[PfP(7r/a)] = +1. Thus, the topological phase 
boundary given by the change in the topological index is 
determined by PfP(O), and, thus, is accompanied by the 
gap closing at zero momentum. Note that the topologi- 
cal reconstruction of the spectrum is always accompanied 
by closing of the bulk gap [9] since Deti^BdC = PfP^. 
Our approach of calculating the TP invariant relies on 
translational symmetry. Recently, the expression for the 
TP invariant was generalized to spatially inhomogeneous 
case [62l [63j. 
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We now calculate FfB{0) for a simplified three- 
subband model and compare the phase diagram with the 
numerical one presented in the previous section. To make 
progress we assume that A^^^ = A, A^^^+i = A' and 



consider only diagonal in the subband index spin-orbit 
coupling terms in Eq. Q. With these approximations, 
FfB{0) becomes 



FfB{0) = {SEf^ {-V^ + A^ 



A'2 + 



^13 



+(5^13 {{V^ - A2)A'2 + (2T/2 _ 2A2 - 3A'2) - 2/."))) 



J -r 2 (A^ + 2A'2 
+ (2A2 + 3A'2 

M^) + M + A^ + (-v;2 + A^ 

.2 



(A- A')^ 
(A' + 



{V^-iA + A') 



2A'' 



m' + /-v;'(2a2 + 3A'2 + 2m')) 



2A'- 



/\2 



(33) 



Here (^£^12(3) represents the energy difference between 
first and second (third) subbands due to the transverse 
confinement. The superconducting gaps are related to 
the nominal bulk gap Aq via relations A = 7Ao/(7 + Aq) 
and A' = 0.257Ao/(7 + Aq) which take into account 
the dependence of the induced parameters on tunnel- 
ing strength. Here we have chosen a reasonable ratio 
of AyA = 0.25. Similarly, all other energies are renor- 
malized in the following way E ^ ^^Aq explained 
in the previous sections. 

The phase diagram showing a sequence of topologi- 
cal phase transitions for the three sub-band toy model is 
shown in Fig. 13 The panels (a)-(c) represent the phase 
diagram with no interband mixing terms (i.e. A' = 0). 
One can clearly see crossings in the phase diagram which 
represent the "sweet spots". One can also notice the ef- 
fect of the renormalization due to SM-SC tunneling as 
we increase 7 - the superconducting and normal terms 
are rescaled in a different way as explained above. In the 
weak-coupling regime 7 <C Aq, the normal terms are not 
significantly renormalized since ^^^^ ~ 1 whereas in- 
duced pairing is small A ^ 7 and is entirely determined 
by the normal state level broadening 7. On the other 
hand, in the strong-coupling regime 7 Aq, the normal 
terms are decreasing with 7 because ^^^o ^ A0/7 <C 1 
whereas A saturates at Aq. Thus, strong SM-SC tunnel- 
ing leads to important quantitative effect which should 
be taken into account in a realistic model for the prox- 
imity effect. 

The panels (d)-(f) represent the phase diagram with 
finite interband mixing terms A' 7^ 0. Here we find qual- 
itative agreement with the numerical results presented 
in the previous sections, compare Figs. |8] and |9] with 



(d) and (e). Different renormalization of the normal and 
SC terms has a two-fold effect on the phase diagram: 

a) with the increase of tunneling the topological phase 
is effectively "pushed" towards higher magnetic fields; 

b) even small interband pairing term opens a large gap 
at the sweet spot leading to extended vertical topolog- 
ical regions. This insensitivity of the topological phase 
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FIG. 13: (Color online) Phase diagram of the multiband 
nanowire at different SM-SC tunneling strengths 7 and in- 
terband pairing A12 = A23 = A^ calculated analytically for 
the effective three-band model, (a)-(c) correspond to 7 = 
5,20,40, respectively, and superconducting gaps A/Ea = 4 
and A'/Ea = 0. (d)-(f) correspond to 7 = 5,20,40, respec- 
tively, and superconducting gaps A/Eo, = 4 and A! IE a — 1. 
We used here 5E12/E0, = 30 and dEis/Ec, = 80. 



against chemical potential fluctuations can be exploited 
for the protection against disorder in the multisubband 
occupied nanowires with no such situation arising in the 
single channel case. 



E. Phase diagram in the presence of a transverse 
external potential. 

A key feature of the phase diagrams shown above is 
represented by the hot spots. In the presence of in- 
ter sub-band pairing, the topologically nontrivial phases 
expand in the vicinity of these hot spots, and, for cer- 
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tain values of the Zeeman field, become stable over a 
wide range of values for the chemical potential. As we 
show in the subsequent sections, this property is criti- 
cal for stabilizing the topological superconducting phase 
and, ultimately, for realizing and observing the Majo- 
rana zero modes. The necessary ingredient for obtaining 
this expansion of the hot spots is a non-vanishing inter- 
sub-band pairing, which is obtained by a non-uniform 
SM-SC coupling characterized by a strong position de- 
pendence along the direction transverse to the wire (see 
Fig. [3|. The natural question is whether this effect 
can be obtained by breaking the symmetry in the trans- 
verse direction using an external field, e.g., generated 
by a gate potential, instead of a non-uniform coupling. 
This would constitute an alternative to engineering non- 
uniform SM-SC interfaces that would be much simpler 
to implement and would allow better control. To inves- 
tigate this possibility, we consider an external potential 
that varies linearly across the nanowire. 



Kxt(y) = 1/o(2y/iy-l), 



(34) 



where Vq is the amplitude of the transverse external po- 
tential. The matrix elements for the external potential 
are strictly off-diagonal and couple strongly the neigh- 
boring sub-bands, while other contributions are at least 
one order of magnitude smaller. Consequently, to a first 
approximation we have , {ny\Ve^t\n'y) ^ 0.41/o(^n;,n^±i. 
Even in the presence of this inter-band coupling, the 
Pfaffian Pf5(0) for the effective three-band model can 
be determined analytically and it is given by an expres- 
sion that generalizes Eq. (33). The corresponding phase 



diagrams for both uniform and non-uniform SM-SC cou- 
plings and different values of the external potential are 
shown in Fig. [T4j 

The key conclusion of this calculation is that inter- 
sub-band mixing due to an external transverse potential 
does not lead to an expansion of the topological phase 
in the vicinity of the sweet spots, but rather determine a 
shift in their position. This is illustrated by diagrams (a- 
c) in Fig. 14, which correspond to homogeneous SM-SC 
coupling, i.e., ^ = 0, and different values of the external 
field. Notice that, to a first approximation, the effect of 
the transverse potential is equivalent to increasing the 
inter-sub-band spacing SEnyu'y- If the transverse poten- 
tial is applied across a wire that is non-uniformly coupled 
to the SC, in addition to shifting the position of the sweet 
spots, it reduces the stability of the topological phase in 
their vicinity, as illustrated in panels (d-f) for ^ 0.4. 
In addition, as a consequence of effectively increasing the 
inter-sub-band spacing, the regions between successive 
sweets spots occupied by the topological SC phase ex- 
pands with increasing the transverse potential, as shown 
in panels (c) and (f). Note that, in the presence of non- 
uniform SM-SC coupling, the phase diagrams for Vb > 
and Vb < (not shown in Fig. 14 are slightly differ- 



ent. The strong dependence of the phase boundaries on 
the transverse potential, especially in the vicinity of the 
sweet spots, could be used experimentally for driving the 
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FIG. 14: (Color online) Phase diagram of a multi-band wire 
in a transverse external field. The strength of the SM-SC 
coupling is 7 = Aq, with uniform coupling {0 — 0, i.e., = 0) 
for the top panels and non-uniform tunneling (^ ^ 0.4, A^ = 
A/4) corresponding to the bottom panels. The amplitude of 
the transverse potential is: Vb = for panels (a) and (d), 
Vb = 50^a for (b) and (e), and Vb = 100^;^ for (c) and (f). 
Note that the topologically nontrivial phase has a vanishing 
width at the sweet spots in the absence of inter-sub-band 
pairing (a-c) for all values of the external potential. For non- 
uniform coupling (d-e), the transverse potential reduces the 
width of the topological phase near the sweet spots. . 



system across the topological phase transition by tuning 
a gate potential, instead of changing the magnetic field. 



IV. EFFECT OF DISORDER ON THE 
TOPOLOGICAL SUPERCONDUCTING PHASE 



In this section we consider the effect of disorder on the 
stability of the topological superconducting phase har- 
boring Majorana fermions. In a realistic system, dis- 
order comes in various ways that affect the topological 
phase very differently. In this paper we consider three 
types of disorder: impurities in the s-wave superconduc- 
tor, short- and long-range disorder in the semiconductor 
wire, and random nonuniform coupling between the semi- 
conductor wire and the superconductor, which mimics a 
rough interface and the imprecision in engineering inho- 
mogeneous ^/-dependent couplings. We begin by consid- 
ering short-range impurities in the bulk superconductor, 
followed by the consideration of the other two types of 
disorder which are both more complex to treat theoreti- 
cally and more detrimental to the existence of the Majo- 
rana modes. 
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FIG. 15: (Color online) Diagrammatic perturbation theory 
in the tunneling between semiconductor and superconductor. 
Disorder-averaging is performed at each order in tunneling 
t. The thick solid line represents disorder- averaged Green's 
function in the superconductor G(jp^uj). The last diagram 
corresponds to irreducible contributions, as far as disorder 
averaging is concerned. 



Short-range impurities in the bulk 
superconductor. 



In order to understand the effect of non-magnetic im- 
purities on the induced superconductivity in the semicon- 
ductor, we first review the results on the proximity effect 
for the infinite planar interface presented in the previous 
section. The basic idea is that the presence of short- 
range non-magnetic impurities in the metal modifies the 
bulk Green's function, which then is used to derive the 
appropriate superconducting proximity effect. 

We begin our discussion by considering the pertur- 
bation theory in the tunneling t, which is justified in 
the limit of low interface transparency. The lowest or- 
der contributions of the diagrammatic expansion in t 

One can notice that the self- 



are shown in Fig. 15 



energy at the second order in t is determined by the 
disorder- averaged Green's function in the superconduc- 
tor. Since typical s-wave superconductors are disordered 
(i.e. tAq <C 1 with r being momentum relaxation time), 
the effect of impurity scattering is non-perturbative. In 
general, this yields a non-trivial problem because of the 
self-consistency condition which now has to be solved 
in the presence of disorder. The problem, however, can 
be substantially simplified if we neglect the effect of the 
in-plane magnetic field on the s-wave superconductivity. 
This condition can be justified due to the vast difference 
of the ^-factors in the superconductor and semiconduc- 
tor. We also assume here that the superconductor is 
thin enough so that we can neglect orbital effects. At 
this level of approximation, the problem reduces to un- 
derstanding the effect of non-magnetic impurities on the 
bulk s-wave superconductivity, which has been investi- 
gated long time ago by Abrikosov and Gor'kov [BJ- The 
main result of Abrikosov- Gor'kov's theory is that dis- 



order does not affect s-wave superconductivity, i.e., the 
superconducting gap is not suppressed by non-magnetic 
impurities, in agreement with Anderson's arguments in- 
voking time-reversal symmetry. Specifically, the disorder 
at the single-particle Green's function level merely leads 
to the renormalization of the parameters Aq = ?7a;Ao and 
Co r]^Ao: 



1 



1 



2t^ 



(35) 
(36) 



Here bar represents disorder-averaged Green's function 
and the disorder strength is parameterized by the impu- 
rity scattering time r defined as 1/r = v{eF)TiiUQ with 
i'{eF)i and uq being the density of states at the Fermi 
level, the impurity concentration, and the scattering po- 
tential, respectively. Thus, at this level of the perturba- 
tion theory, the proximity effect can be included using the 
formalism developed for the clean case, see Eq. (13). By 



integrating out the degrees of freedom corresponding to 
the superconductor, one obtains the interface self-energy: 



S(cj) = \t\'u{eF) / deG{e,uj 



= -\t\MeF) 



(37) 



(38) 



Finally, upon substituting the expressions for Co and Aq, 
we find that proximity-induced superconductivity is not 
affected by disorder in the s-wave superconductor 



UJTq + AqTx 



C 



(39) 



Moreover, one can see that, unlike impurities in the semi- 
conductor, the disorder in the superconductor does not 
lead to momentum relaxation in the semiconductor at 
this order of the perturbation theory. 

As the interface transparency is increased, one needs 
to consider higher-order terms in tunneling. These 
terms involve reducible and irreducible contributions, 
see Fig. 15 The former depend only on the disorder- 
averaged Green's function, and, in a sense, are easy 
to take into account, whereas the latter involve non- 
trivial higher-order correlation functions (i.e. diffusons 
and Cooperons) and lead to momentum relaxation in 
the semiconductor . In this work, we consider only re- 
ducible contributions and neglect higher order correlation 
functions with respect to disorder- averaging which, one 
can show, are much smaller than the reducible ones [65j. 
Therefore, our minimal treatment of the disorder in the 
superconductor is justified and we believe that the su- 
perconducting disorder is irrelevant for the topological 
superconductivity. However, the disorder in the semi- 
conductor (and at the interface) is relevant, as we discuss 
next. 
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This conclusion holds as long as we neglect the effect of 
external magnetic field on the s-wave superconductivity, 
which is allowed as long as the field is not too large to 
destroy the superconducting state. This is realistic due 
to the large difference between the ^-factors in supercon- 
ductors (^sc ~ 1) and semiconductors (^sm 10 — 70), 
allowing one to always find a parameter regime where 
the magnetic field opens a large spin gap in the spec- 
trum without destroying superconductivity. Also, note 
that our conclusion regarding the short-range impurities 
in the superconductor is valid when the strength of the 
disorder is much larger than the superconducting gap: 
tAq ^ 1 but is small enough not to affect the density 
of states. We have restricted our analysis to the experi- 
mentally relevant regime 7 < Aq. In the opposite limit 
7 ^ Ao, electrons spend most of the time in the s-wave 
superconductor where their dynamics is not governed by 
the helical Hamiltonian required to have a topological 
phase and, thus, this limit is not experimentally desir- 
able, see discussion in Sec. II B 1| 

We emphasize that our finding that short-range impu- 
rities in the bulk superconductor do not affect the topo- 
logical superconducting phase emerging in the semicon- 
ductor is quite general since the electrons in the nanowire, 
being spatially separated from the bulk superconductor, 
simply do not interact directly with the short-range dis- 
order in the bulk superconductor. Our use of the short- 
range impurity model to characterize the disorder in the 
bulk superconductor is justified, since strong metallic 
screening inside the bulk superconductor would render 
all bare long-range disorder into screened short-range 
one. Thus, as long as the applied magnetic field does 
not adversely affect the s-wave superconductivity, our 
conclusion regarding the validity of the Anderson the- 
orem (i.e., no adverse effect from non-magnetic impu- 
rities) to the whole superconductor-semiconductor het- 
erostructure system applies. In this context, we men- 
tion the recent theoretical analysis of Ref. [66 , where 
it was explicitly established that, in the structure we are 
considering, the proximity-induced superconducting pair 
potential remains s-wave both in the topological and in 
the non-topological phase, in spite of the non- vanishing 
Zeeman field. The existence of the s-wave pairing po- 
tential, even in the presence of a parallel magnetic field 
(provided it is not too large), is the key reason for the 
short-range disorder in the superconductor not having 
any effect on the topological phase. 



B. Disorder in the semiconductor nanowire 

In sharp contrast to disorder in the superconductor, 
disorder in the nanowire can have significant adverse ef- 
fects on the stability of the topological SC phase. There 
are several different potential sources of disorder in the 
semiconductor. We focus on two sources that are the 
most relevant experimentally: random variations of the 
width of the SM nanowire and random potentials cre- 



ated by charged impurities. The first type of disorder 
is generally long-ranged, while the second type can be 
either long- or short-ranged. How to account for the 
effects of disorder depends crucially on the type of phys- 
ical quantity that one is interested in. Here we focus 
on the low-energy spectrum, which carries information 
about the stability of the Majorana bound states, and 
on thermodynamic quantities such as the local density 
of states, which are experimentally relevant. The effect 
of disorder on these quantities is sample-dependent. We 
emphasize that averaging over different disorder realiza- 
tions is equivalent in this case with sample averaging. As 
the goal is to observe stable Majorana fermions in a given 
nanowire, we investigate here the spectrum of the system 
for a given disorder realization and focus on establish- 
ing the general parameter regimes (e.g., amplitudes and 
length scales of the disorder potential) consistent with re- 
alistic experimental conditions that ensure the stability 
of the topological SC phase. More specifically, we study 
several different disorder realizations characterized by a 
given set of parameters and extract the generic features 
associated with that type of disorder. 

For a single channel Majorana nanowire, one can ob- 
tain some analytical results by performing disorder av- 
eraging [67H69] . Specifically, for a model of the spin- 
less p-wave superconductor it has been shown [67l [68] 
that disorder drives the transition into non-topological 
phase when impurity scattering rate becomes compara- 
ble with the induced superconducting gap. In more real- 
istic spinful models involving semiconductor nanowires, 
the physics is richer and depends on the strength of the 
magnetic field. We refer the reader to Ref. [69 for more 
details. The generalization of these results to the case 
of disordered multi-band semiconductor nanowires is an 
interesting open problem. 



Semiconductor nanowires with random 



The dimensions of the nanowire in the transverse direc- 
tion satisfy the relation Ly ^ Lz. The small thickness 
Lz is critical for the effectiveness of the superconduct- 
ing proximity effect, as the SM-SC effective coupling 7 
scales, approximately, as l/L^z- On the other hand, the 
much larger width Ly is required by the multichannel 
condition. Atomic-scale variations of Lz generate huge 
local potential variations (of the order 500£^q,) that would 
effectively cut the wire in several disconnected pieces. 
Topological SC phases may exist inside each of these 
pieces, but the Majorana states will be localized at the 
boundaries separating different segments and, in general, 
tunneling between them will be nonzero. To realize a sin- 
gle pair of Majorana zero-energy states localized at the 
ends of the wire, Lz should be uniform along the system. 
(We mention in passing that modern MBE growth is con- 
sistent with very small variations in Lz as necessary for 
the realization for the Majorana. ) 

Engineering a long wire with constant width Ly may 
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FIG. 16: (Color online) Broadening of the sub-bands in a 
nanowire with random edges. For Uy > 2 the broadening 
becomes comparable with the inter-band gap in the presence 
of fluctuations AL^ representing a few percent of the width 
Ly, i.e., the sub-bands loose their identity. 



be, on the other hand, extremely challenging and less 
relevant for the stability of the topological SC phase. 
We assume that Ly varies along the wire randomly in 
atomic-size steps that extend along hundreds of lattice 
sites in the x-direction, resulting in a function Ly{x) that 
varies over length scales larger than the width of the wire. 
These fluctuations generate local variations of the bare 
sub-band energies given by Eq. ([3| of the order of 
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where Ly is the average width of the nanowire, AL^y the 
value of the local variation, and a the lattice spacing. 
Note that different Uy sub-bands are shifted differently, 
i.e., the random edge is not equivalent to long-range 
chemical potential fluctuations. To provide a quantita- 
tive measure of this effect, we show in Fig. [16] the evolu- 
tion of the nanowire sub-bands with the size of the fluctu- 
ations /S.Ly. Note that the sub-bands loose their identity 
in the presence of fluctuations representing a few percent 
of the wire width, as the broadening becomes compa- 
rable with inter-sub-band gaps. The natural question 
is how this broadening affects the low-energy physics of 
the nanowire and, in particular, the topological SC phase. 
Intuition based on the weak-coupling picture would sug- 
gest that the topological phase might become unstable, 
as the parity of the number of sub-bands crossing the 
chemical potential becomes an ill defined quantity. 

To quantify the effect of random edges on the low- 
energy physics, we first parametrize this type of disorder. 
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FIG. 17: (Color online) Profile of the variation of the nanowire 
width, ALy, as a function of position along the wire. The 
characteristic length scale for profile I (upper panel) is about 
10% of the nanowire length L^, wile for profile II (lower panel) 
it is approximately 5%. Each profile generates a series of par- 
ticular disorder realizations characterized by different values 
of the maximum amplitude |ALy|max- Random edge pro- 
files corresponding to a given characteristic length and hav- 
ing the same maximum amplitude have similar effect on the 
low-energy physics of the nanowire. 



We consider a nanowire of width Ly{x) = Ly -\- AL^(x), 
where ALy{x) is a random function characterized by by 
a certain maximum amplitude and a characteristic wave- 
length. Two examples of random edge profiles are shown 
in Fig. 
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We note that the actual width of the wire 
varies in atomic steps, i.e, lAL^ylmin = cl and we assume 
that the characteristic length scale of these variations is 
much larger than the atomic scale. For example, for a 
nanowire with L^ = 5/im, a random edge profile like in 



the upper panel of Fig. 17 and a maximum amplitude of 
10%, the atomic edge steps extend over hundreds of unit 
cells. In the calculations we explore the effect of random 
edges with maximum amplitudes up to 10% of Ly and 
various characteristic wavelengths. 

The mini-gap A* is reduced by the presence of random 
edges. However, for control parameters corresponding to 
points in the phase diagram away from phase boundaries, 
the amplitude of the variations of Ly required for a com- 
plete collapse of the gap is well above 10%. Examples of 
the spectra for systems with random edges are shown in 
figures 18 and 19 Based on a number of similar calcu- 



lations for different disorder realizations and control pa- 
rameters, we have established the following general con- 
clusions: 
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FIG. 18: (Color online) Spectrum of a weakly coupled 
nanowire (7 = 0.25Ao) with random edges in the topological 
phase with = 1 (see Fig. [8|. In the presence of disorder the 
mini-gap A* decreases, but remains finite. The calculation 
was done using the random edge profiles shown in Fig. |17| for 
maximum amplitudes of 2% and 10%. Longer range disorder 
(upper panel) has stronger effects than short range disorder 
(lower panel). Small amplitude variations of the nanowire 
width of the order 2% generate a reduction of A* up to 30%, 
but further increasing the amplitude has a weak effect at low 
energies. 



i) The details of the low-energy spectrum of a dis- 
ordered nanowire depend on the particular disorder re- 
alization. Nonetheless, different disorder profiles charac- 
terized by a given amplitude and having the same charac- 
teristic length scale are likely to generate similar values 
of the mini-gap A*, with the exception of a few "rare 
events", which are characterized by significantly lower 
gap values. A calculation that involves averaging over 
disorder will capture these "rare events" and will predict 
a value of the gap much lower than the typical value. 
Such a calculation would be relevant for an extremely 
long wire, i.e., in the limit ^ 00 , or for systems with 
a very large density of states at the relevant energies (e.g., 
a metal). However, in a typical semiconductor wire the 
number of states that control the low-energy physics is 
of the order of 100. How the energies of these states 
are modified in the presence of disorder, depends on the 
specific details of the disorder profile. Hence, any exper- 
imentally relevant conclusion regarding the low-energy 
spectrum or the local density of states of a disordered 
nanowire should be based on calculations involving spe- 
cific disorder realizations. A direct consequence of these 
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FIG. 19: (Color online) Spectrum of a nanowire with random 
edges at intermediate coupling, 7 = Ao . The control param- 
eters correspond to a point in the phase diagram inside the 
topological phase with N = 1 (see Fig. |9|. The random edges 
are given by the profiles in Fig. fTTl 



considerations is that nominally identical samples with 
the same average disorder (e.g. same mobility) may have 
very different Majorana minigaps since they are likely to 
have different disorder configurations. The distribution 
of the minigaps was recently studied in Ref. [70 . 

ii) Long range disorder has a stronger effect than short 
range disorder. This is general characteristic of disor- 
dered nanowires, regardless of the source of disorder, and 
will be studied in more detail in the next section. We 
note that the effects of variations of wire width ALy{x) 
at atomic length scales are negligible. 

iii) Intermediate coupling 7 ~ Aq represents the opti- 
mal SM-SC coupling regime. The large gap that char- 
acterizes the topological = 1 phase in this regime is 
robust against fluctuations of the wire width of the order 
±10% for any set of parameters that are not in the im- 
mediate vicinity of a phase boundary. Most importantly, 
this condition can be satisfied for specific values of the 
Zeeman field (e.g., of the order SOE^ for the parameters 
corresponding to the phase diagram in Fig. |9| over a 
large range of chemical potentials. 



2. Nanowires with charged impurities 

A major source of disorder in the nanowire consists of 
charged impurities. Because the carrier density in the SM 
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nanowire is small, charged impurities are not effectively 
screened and their presence can potentially have signifi- 
cant effects. The nature of these charged impurities, the 
values of their effective charge and their exact locations 
depend on the details of the nanowire engineering process 
and will have to be determined by future experimental 
studies. Here, we are interested in the fundamental ques- 
tion regarding the stability of the topological SC phase. 
In particular we address the following question: Is it pos- 
sible to realize stable zero-energy Majorana modes in a 
nanowire with charged impurities within a realistic sce- 
nario? To answer this general question, we focus on four 
key aspects of the problem: a) The screening of charged 
impurities by the electron gas in the SM nanowire, b) The 
dependence of the low-energy physics on the concentra- 
tion of impurities, c) The dependence of the low-energy 
spectrum on the Zeeman field, and d) The effect of long 
range random potentials. 

a. Screening of charged impurities. We start by con- 
sidering a single charge q inside or in the vicinity of the 
nanowire. For concreteness we assume that the charge 
is positioned near the middle of the wire and one lat- 
tice spacing away from its surface, i.e., for a wire that 

< Xn < Ln, with 



the position of the impurity is given by 
zmp) = {Lxf^^ —a). This corresponds. 



occupies the volume defined by 
j e {x,y,z} 

for example, to a charged impurity localized at the in- 
terface between the SM and the SC. We consider the 
extreme case g = ±e, where e is the elementary charge, 
although in practice it is likely that the effective charge 
is only a fraction of this value due to screening by elec- 
trons in the SC. We neglect the screening due to the 
presence of the superconductor, which may significantly 
reduce the effective potential created by the charge. A 
simple estimate of the screening effects due to the elec- 
trons in the semiconductor within the Thomas-Fermi ap- 
proximation is highly inaccurate due to the low carrier 
density. We checked this property explicitly by calculat- 
ing numerically the carrier density induced by a given 
effective potential V{r), e.g., a screened Coulomb po- 
tential. We find that the relationship between the in- 
duced local carrier density Sn{r) and the local effective 
potential is highly non-linear. In addition, the density 
is characterized by strong Friedel-type oscillations (see 
Fig. 20). Hence, solving quantitatively the screening 
problem for the nanowire would require a self-consistent 
calculation that includes electron-electron interactions at 
the Hartree-Fock level. This calculation is beyond the 
scope of the present study and will be addressed else- 
where. Here we address a more limited question: What 
is the characteristic length scale over which the exter- 
nal charge q is screened? We define this length scale A as 
the characteristic length of the volume that contains 63% 
(i.e., a fraction equal to 1 — 1/e) of the induced charge. 
Specifically, the potential created by the charge at a 
point r inside the semiconductor is 
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FIG. 20: (Color online) Induced carrier density as a function 
of distance from the impurity. The two curves, corresponding 
to different external sets of parameters for the nanowire, have 
been shifted for clarity. The upper curve corresponds to a 
system with a single occupied sub-band and is characterized 
by A 45a (i.e., 63% of the induced charge is in a disk of 
radius 45a), while the lower curve is for a system with three 
occupied sub-bands and has A ^ 30a. Note that the effective 
impurity potentials that generate these density profiles are 
given by Eq. (42) with the corresponding values of A. 



where rimp = (^imp, ^imp, ^imp) is the position of the 
impurity and is the relative dielectric constant of the 
SM. For In As = 14.6. Next, we take into account the 
fact that the nanowire is extremely thin in the z direction 
{Lz ~ 10a), hence the wave function profile along this 
direction is very little affected by the presence of the 
impurity. On the other hand, the induced charge has a 
strong dependence on x (the direction along the wire) 
and y (the transverse direction). As mentioned above, 
the impurity potential is screened by the induced charge 
outside a region with a characteristic length scale A that 
contains a fraction equal to 1 — 1/e of the induced charge. 
We assume that the screened potential is qualitatively 
described by an expression of the form 



Vs = V{r - rimp)exp 



y^{x- Ximp^ ^{y- Vimpy 



A 



(42) 

The parameter A from Eq. (42) is determined self- 
consistently by imposing the condition that 63% of the 
induced charge be in a disk of radius A and thickness 
Lz centered at {ximp^Vimp^ Lz/2). The results for two 
different sets of parameters are shown in Fig. 
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(41) 



expected, nanowires with multiple occupied sub-bands 
provide a more effective screening, which is reflected in 
a lower value of the screening length A. We empha- 
size that the present approach is not fully self-consistent 
and is therefore only of qualitative validity. An effective 
screened potential that includes exactly the contribution 
of the induced charge with details (e.g., oscillatory com- 
ponents) that are not captured by Eq. ( [42| should give 
results very similar to what we obtain here. Nonethe- 
less, we expect these details to have a weak effect on 
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FIG. 21: (Color online) Top panel: Comparison between the 
low-energy spectrum of a clean nanowire (yellow squares) and 
the spectrum of a nanowire with a charged impurity with q = 
e (red circles). Note that the two spectra are almost on top of 
each other. Bottom panel: Electron wave function amplitudes 
for two low-energy states (marked by arrows in the top panel) 
in the presence of the charged impurity. The Major ana zero 
modes (n = ±1, no t shown) and the in-gap modes (n 
±2, ±3, see Fig. 
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and the corresponding discussion) are 
localized near the ends of the wire and are not affected by the 
presence of the impurity. The low-energy bulk states |n| > 4 
that are extended in a clean nanowire become localized near 
the impurity. 



the final results. We also note that the self-consistent 
calculation of the screening length is done for a non- 
superconducting nanowire, then the effective impurity 
potential Vg is added to the total nanowire Hamiltonian 
before including the proximity effects due to the SM-SC 
coupling. 

What is the effect of the charged impurity on the low- 
energy spectrum of the superconducting nanowire? The 
effective potential in the vicinity of the impurity is ex- 
tremely large, e.g., Vs{Lx/2^Ly/2^Lz/2) ^ 640meV (i.e., 
^ 12000£^Q,), much larger than any other relevant energy 
scale in the problem, hence one would naively expect sig- 
nificant effects. However, the low-energy physics of the 
SM nanowire is controlled by single particle states with 
low wave numbers and the matrix elements of the impu- 
rity potential between these state are relatively small. In 
other words, the impurity will strongly affect the spec- 
trum at intermediate and high energies, but will have a 
relatively small effect at low energies. To illustrate this 
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FIG. 22: (Color online) Specific disorder realizations in 
nanowires with charged impurities. The dots represent the lo- 
cations of the impurities. Each impurity has a charge q — =be 
and is positioned one lattice spacing away from the surface 
of the nanowire along the z direction. The linear impurity 
densities riimp are: l/j.m~^ (A), 2/j.m~^ (B), 4/j.m~^ (C), and 
(D). 



low-energy spectra of a clean nanowire and of a nanowire 
with a charged impurity. We note that the presence of 
the external charge induces the formation of localized 
states in the vicinity of the impurity (see Fig. 21). 



property, we show in Fig. 21 a comparison between the 



Charge impurities in one-dimensional quantum wires 
produce weakly long-range disorder since the Coulomb 
potential decays as Inlg'al in the momentum space for 
q ^ with a being the short-distance cut-off associated 
with the transverse dimensions of the nanowire [71j . This 
is to be contrasted with the much stronger q~^ {q~'^) ^ong 
wavelength divergence of the bare Coulomb disorder in 
two (three) dimensions. The weakly long-range nature 
of ID Coulomb potential suggests that any regulariza- 
tion of the long-range disorder would be a reasonable 
approximation in spite of the fact that Thomas-Fermi 
screening itself is weak in ID. In particular, the pres- 
ence of interband scattering in the multisubband situa- 
tion would essentially lead to effective 2D screening in 
the system, which should suffice to regularize the sin- 
gular Coulomb disorder. At very low densities, where 
the nanowire is strictly in ID limit with only the low- 
est subband occupied, weak screening will lead to the 
formation of the inhomogeneous electron puddles in the 
system around the charge impurities due to the failure of 
screening. This situation is detrimental to the Majorana 
formation and must be avoided. It is clear that higher 
density and multisubband occupancy would be favorable 
for the experimental realization of the Majorana modes 
in the semiconductor nanowires. 

h. Charged impurity disorder. The next question 
that we address concerns the dependence of the low en- 
ergy spectrum of a nanowire with charged disorder on 
the concentration of impurities. We emphasize that the 
details of the low-energy physics depend on the specific 
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disorder realization. In particular, the low-energy states 
are localized in the vicinity of the impurities (see Fig. 
21 ) and their energies depend on the specifics of the real- 
space disorder configuration. Hence, as discussed above, 
calculations of single-particle quantities, e.g., the local 
density of states, should involve specific disorder realiza- 
tions, rather than disorder averaging. In Fig. [22]we show 
four different specific disorder realizations, correspond- 
ing to linear impurity densities ranging from l/im~^ to 
8/im~^, which are reasonably realistic impurity densities 
(~ 10^^ cm~^) in high quality semiconductor structures. 
The impurities carry charge q = e and are positioned 
at a distance of one lattice constant away from the SM 
surface. The effect of these impurities is incorporated 
through an impurity potential of the form 



Vimp{r) = ^V.,{r,rj), 



(43) 



where Vs is given by Eq. (42) and Vj are the impurity 



position vectors. We note that in Eq. (43) the screened 



potential is characterized by a screening length A deter- 
mined as described above for a single impurity. This ap- 
proximation does not take into account the impact of the 
dependence of the effective potential associated with a 
given charged impurity on the location of the charge and 
on the presence of other impurities. It also neglects the 
effect of screening by the SC itself which should strongly 
suppress the effective disorder. 

The low-energy spectra of a nanowire with random 
charged impurities distributed as in Fig. [22] are shown in 
Fig. 23 for two sets of control parameters. The general 
trend is that the mini-gap decreases with increasing im- 
purity concentration. However, for a given concentration 
nimp the exact value of the mini-gap depends on the spe- 
cific disorder realization. As mentioned above, averaging 
over disorder includes rare configurations characterized 
by small mini-gaps, hence the averaged density of states 
is characterized at low energies by a small weight that 
does not correspond to any physical state in a typical 
disorder realization. The signature of a topological SC 
phase with N = 1 (i.e., one pair of Majorana fermions) 
that distinguishes it from the trivial SC phase with N = 
is the presence of zero-energy quasiparticles separated by 
a finite gap from all other excitations. One key conclu- 
sion of our calculations is that the mini-gap that protects 
the topological SC phase remains finite for a significant 
range of realistic impurity concentrations. 

c. Dependence of the low-energy spectrum on the 
Zeeman field. As shown in Fig. [23j for certain high im- 
purity concentrations (e.g., configuration D in fig. 22) it 
is possible that, in addition to the Majorana zero mode, 
another low-energy state has energy close to zero and is 
separated by a gap from the rest of the spectrum. How 
can one distinguish experimentally this state from a topo- 
logical superconductor characterized by a finite mini- 
gap, on one hand, and from a trivial superconductor with 

= 2, on the other? As the low-energy spectrum has a 
strong dependence on F, the key is to vary the Zeeman 
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FIG. 23: (Color online) Low-energy spectra of a supercon- 
ducting nanowire with charged disorder for two different sets 
of control parameters and four disorder realizations. The sym- 
bols correspond to the specific disorder realizations shown 
in Fig. 22 Note that: i) The mini-gap A* finite for 
impurity concentrations nimp < 4/im~^ and collapses for 
^imp = 8/i?n~^, and ii) There is an overall tendency of the 
low-energy features to move at lower energies when increas- 
ing the impurity concentration, but the exact value of the 
mini-gap depends on the specific disorder realization, e.g., in 
the lower panel A*(nimp = 4) > A*(nimp = 2). 



field. Tuning F may push the system into the topologi- 
cal SC phase and the energy of the additional state will 
increase. On the contrary, it is possible that varying the 
Zeeman field will lead to the appearance of more low- 
energy excitations. This is the characteristic signature of 
the transition zone between phases with different topolo- 
gies. To address this problem more systematically, let 
us consider a nanowire with the same parameters as in 



the bottom panel of Fig. [Toj fi = 54.5£^q,, F = Aq, and 
6 = 0.8. The dependence of the quasiparticle gap on F 
for the infinite clean system in shown in Fig. [lOj and 
the dependence of the mini-gap on F for a finite wire 



is shown in Fig. 12 The transition between the trivial 
SC phase with = and the topological phase with 
N = 1 is clearly marked by the vanishing of the gap at 
F ^ 21.bEa. 

Let us now add disorder and follow the evolution of 
the lowest three energy levels with the Zeeman field. The 
results for two different values of the impurity concentra- 



23 




10 15 20 25 30 35 




10 15 20 25 30 35 

Zeeman field F/Ea 

FIG. 24: (Color online) Dependence of the lowest three en- 
ergies En {n = 1,2,3) on the Zeeman field for a disordered 
nanowire with riimp — Aiim~^ (top) and riimp — 7/j.m~^ (bot- 
tom). The lowest energy state is characterized by a gap = 
El (red/dark gray region) that vanishes for F > 21.6^a, i-e., 
when the system enters the topological SC phase. The gap be- 
tween the first and the second levels, A 2 = E2 — Ei (yellow), 
becomes the mini-gap in the topological phase. The tran- 
sition zone is characterized by a high density of low-energy 
modes and expands with increasing the impurity concentra- 
tion. 



tion are shown in Fig. 24 First, we note that the main 
features are similar with those observed in a clean system: 
at low values of the Zeeman field the system is in a triv- 
ial SC phase characterized by finite gap to all excitations, 
while for F above a certain critical value Ec ~ 22£^ck the 
system has a Majorana zero mode separated by a mini- 
gap from the rest of the spectrum. The major difference 
from the clean case consists of the transition zone, which 
extends over a finite range of values of the Zeeman field 
and is characterized by multiple low-energy excitations. 
This transition zone extends with increasing impurity 
concentration. Assuming that we probe the low-energy 
properties of the system with a certain finite resolution, 
e.g., AE = 0.2£^Q,, the topological phase can be unam- 
biguously distinguished from the trivial SC provided the 
mini-gap is larger that the energy resolution. The trivial 
N = phase will be characterized by a finite gap and no 
zero-energy excitation, while the topological phase will 
have the characteristic zero-mode separated from the fi- 
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FIG. 25: (Color online) Low-energy spectra of a nanowire 
with long-range disorder. The system in characterized by 
jj = 54:.5Ea and 7 = Ao and the impurity potential is given 
by Eq. (44 ) with with a profile ^;imp as in Fig. 17 (lower panel) 
and different amplitudes Uq. The parameters Uo and F are 
expressed in units of Ea. For F = 29 the mini-gap collapses 
for Uo > 60, while for a Zeeman field F32 it survives up to an 
amplitude Uo ~ 100. 



nite energy excitations by the mini-gap. In between the 
two phases there will be a transition zone characterized, 
within our finite energy resolution, by a continuum spec- 
trum. Starting with low values of the Zeeman field, i.e., 
deep inside the N = phase, by increasing F one first 
reaches the transition zone, the topological N = 1 phase. 
We emphasize that for A* < AE the topological phase 
becomes indistinguishable from the transition zone. In 
addition, at large impurity concentrations the mini-gap 
will collapse completely. 

d. Long range disorder potential What is the effect 
of a long-range disorder potential on the stability of the 
topological SC phase? Are there qualitative differences 
from the short-range case discussed above? We are not 
interested here in the possible source of such long range 
disorder, but rather in identifying the magnitude of the 
amplitude of the random potential that would destroy 
the topological phase. 

Let us consider a random potential with a character- 
istic length scale X > Ly. Neglecting the dependence of 
the potential on y and 2:, we have 



Vimp{r) = UoVimp{x/La,), 



(44) 



where Uq is the amplitude of the potential and —1 < 
^imp < 1 is a random profile. For concreteness we con- 



sider the profile shown in the lower panel of Fig. 17 and 
a nanowire with = 5 jam and Ly = 0.12 jam. These 
parameters correspond to A ~ 2Ly. The correspond- 
ing low-energy spectra for different values of the disor- 
der amplitude are shown in Fig. 25 A striking feature 
is represented by the significant difference between the 
critical disorder amplitudes at which the mini-gap col- 
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lapses for the two slightly different values of the Zeeman 
field. Note that both values are near the "center" of the 
topological = 1 phase, as one can see, for example, by 
inspecting the bottom panel of Fig. [12] What sets the 
scale for the critical amplitude? A hint can be obtained 
from the phase diagram shown in Fig. [9] In essence, 
in the presence of a smoothly varying random potential 
topological superconductivity is stable as long as the lo- 
cal chemical potential at any point within the nanowire 
has values within the topological phase. The permis- 
sible range of variation for the local chemical potential 
depends strongly on the Zeeman field. From the phase 
diagram in Fig. |9]and the chemical potential dependence 
in Fig. 11 it is clear that F ^ 32Ea represents a value of 
the Zeeman field that allows for large amplitude chemical 
potential fiuctuations. 

We conclude that long range disorder can be treated as 
local chemical potential fiuctuations. We emphasize that 
this is not the case for short range disorder. The topologi- 
cal phase is stable against chemical potential fiuctuations 
up to a maximum amplitude that depends on the Zeeman 
field, on the average chemical potential, and on the SM- 
SC coupling. Single band occupancy necessarily limits 
the critical amplitude due to the close proximity of the 
phase boundary. Multi-band systems avoid this problem 
and provide a more efficient screening for the short range 
potentials created by charged impurities. In addition, a 
strongly non-uniform coupling between the nanowire and 
the s-wave SC (with ^ ~ 1) with a coupling strength in 
the intermediate regime (7 ~ Aq) provides the optimal 
shape of the phase diagram (see Figures [9] and 12). 



Disorder at the semiconductor— superconductor 
interface 



The inhomogeneous random coupling at the 
semiconductor-superconductor interface is another 
significant source of disorder. In principle, the tunneling 
matrix elements t(ix-,iy) between the nanowire and 
the s-wave SC are characterized by random real space 
variations due to inhomogeneities in the tunneling 
barrier. For example, realizing a nonuniform coupling 
t{iy)^ which is critical for generating off-diagonal pairing 
and for stabilizing the topological phase near the sweet 
spots, may require the growth at the interface of an 
insulating layer with variable thickness across the wire. 
Any growth imperfection will translate into variations 
of t. While ultimately the details of these variations will 
have to be determined by a careful experimental study 
of the interface, it is reasonable to assume that a typical 
interface is characterized by atomic size variations with 
a characteristic length scale of a few lattice spacing, as 
well as longer range inhomogeneities with characteristic 
length scales comparable to the width of the wire, Ly^ 
or larger. The short range inhomogeneities could be 
generated by impurities or by point defects present at 
the interface, while longer range inhomogeneities could 



be due to extended defects. In the absence of a detailed 
microscopic description of the SM-SC interface, it is 
difficult to estimate the amplitude of these ffuctuations. 
Here we consider a phenomenological model of the 
interface and we assume that the tunneling matrix is 
given by 



(45) 



where t{iy) is the smooth component of the nonuniform 
coupling, e.g., the interface transparency with the pro- 
file shown in Fig. [sj and At{ix,iy) represents the ran- 
dom component. To model the short range disorder, we 
coarse grain the interface in square patches of side length 
I and assume that At is uniform within a patch, but 
varies randomly from patch to patch with an amplitude 

(At)max, i.e., with -(At)max < At{ix,iy) < (At)max. 

In the numerical calculations we considered patches of 
sizes / = 20a and I = 60a ~ Ly/A and an amplitude 
(At)max = 0.25t(0), where t{0) is the maximum value of 
the smooth nonuniform tunneling component shown in 
Fig. [3) We note that these are extremely large fiuctu- 
ations of t(ix-,iy)^ larger that the minimum value of the 
smooth tunneling component, t{Ly) = 0.2t(0), which re- 
sult in significant variations of the local effective coupling 
^{ix^iy) cx \t{ix,iy)\'^ . Examples of short-range random 



couplings within the patch model are shown in Fig. 26 



Before we present the results of the numerical calcula- 
tions, we would like to emphasize the specific way that 
interface disorder enters the effective Hamiltonian. While 
the effect of charged impurities can be included through a 
random potential, variations in the SM-SC coupling gen- 
erate randomness in the effective SC order parameter. 
An n' 7 as well as ffuctuations of the renormalization ma- 



trix Z^/^. From equations (25) and (27) we notice that 



the short range ffuctuations of t and, implicitly, the short 
range fiuctuations of 7, are significantly reduced when 
taking the matrix elements with the eigenstates \ny). As 
we are mainly interested in systems with only a few oc- 
cupied sub-bands, integration over y effectively averages 
out fluctuations with characteristic length scales I <^ Ly. 
As an illustration of this property, we consider the case 
of random coupling with / = 60a and (At)max = 0.25t(0) 
for a nanowire with /i = 54. 5^^, F = 32Ea, and 7 = Aq. 
In spite of the relatively large length scale of the fluctua- 
tion, I = Ly/4:^ the variations of the induced gap Aii{x) 
along the wire are only of the order of 10% of the average 
value. The dependence of the induced gap An on the po- 
sition along the wire is shown in the lower panel of Fig. 



20 The amplitude of the An fluctuations is further re- 
duced if we consider shorter range coupling fluctuations. 
A similar behavior characterizes the renormalization ma- 
trix Z^/^. In addition, as the low-energy properties of 
the system are determined by single particle states with 
small wave vectors kx, we expect a further reduction of 
the effect of short range ffuctuations as a result of the 
integration over x. In particular, if the clean, inflnite 
wire is characterized by a maximum Fermi wave vector 
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FIG. 26: (Color online) Random coupling at the SM-SC in- 
terface. The coupling t contains a random component At 
that is constant within patches of length /, but takes a ran- 
dom value within each patch. The patch sizes are / = 20a (top 
panel) and / = 60a (middle panel). The strength of At within 
the patches is color coded. Note that only part of the inter- 
face is shown, as the typical length of a nanowire used in the 
calculations is Of the order lO^a. Bottom: Dependence of the 
induced gap An on the position along the wire for F = 32Ea, 
jjL — 54.5£^a, 7 — Ao, — 0.8, and a random coupling with 
/ = 60a and an amplitude (At)max = 0.25t. The huge local 
variations of the coupling are strongly reduced by the integra- 
tion over y (see main text) and generate fluctuations of the 
order of 10% in Aii(x). 



kp^we expect the low-energy physics to be insensitive to 
random variations of the SM-SC coupling with charac- 
teristic length-scales / < l/kp. 

The effects of a random SM-SC coupling on the low- 
energy spectrum of the nanowire are illustrated in Fig. 
[27| Remarkably, short range fluctuations with ampli- 
tudes up to 25% of the average coupling ai y = (i.e., the 
maximum value of the coupling in a nonuniform profile - 
see Fig. [3| do not destroy the topological SC phase. The 
relatively weak effect of these strong fluctuations is due 
to the implicit averaging involved in the calculation of 
the matrix elements between single particle states with 
low wave vectors. Increasing the characteristic length 
/ makes this type of disorder more effective, as evident 



from the upper panel of Fig. Hence the natural ques- 
tion: what is the effect of long-range SM-SC coupling 
fluctuations? We consider a smooth variation of t along 
the wire with a profile as shown in the upper panel of Fig. 
pT} An amplitude of these fluctuations equal to 25% of 
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t{y = 0) results in the collapse of the gap (see Fig. 
lower panel). However, the topological phase is robust 
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FIG. 27: (Color online) Spectra of a nanowire with random 
SM-SC coupling. Upper panel: Nanowire with short range 
fluctuations of the SM-SC coupling. The random coupling 
is considered within the patch model described in the text 
and corresponds to the distributions shown in Fig. [26] Lower 
panel: Nanowire with long-range SM-SC coupling fluctua- 
tions. The variations of t along the wire have a profile as 

and an amplitude At. 
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shown in the upper panel of Fig 
Note that for At = 0.25t{y = 0) the gap collapses. 



against long range coupling fluctuations with amplitudes 
smaller than 20%. We note that long range variations of 
the coupling strength could result from the engineering 
process of the nonuniform interface. Limiting the ampli- 
tude of these fluctuations below a certain limit of about 
10 — 15% of the maximum coupling strength should be 
a priority of the experimental effort for realizing a topo- 
logical SC state using semiconductor nanowires. 



V. EXPERIMENTAL SIGNATURES OF 
MAJORANA BOUND STATES 

Probing unambiguously the presence of Major ana 
bound states in the superconducting nanowire represents 
a critical task. In this section we show that local spec- 
tral measurements provide a simple and effective tool for 
accomplishing this task. We focus on the local density of 
states (LDOS), which could be measured using, for ex- 
ample, scanning tunneling spectroscopy (STS), and on 
the differential conductance associated with tunneling 
into the ends of the wire. We establish that these mea- 
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FIG. 28: (Color online) Typical spectra for a trivial SC with 
N = (red squares), a topological SC with = 1 (orange 
diamonds), and a trivial SC with N — 2 (black circles). The 
system is characterized by /x = 54. 5^^, = 0.8, and 7 = 
Aq. The finite energy in-gap states (for Gamma — 15^a 
and r = 45^a) together with the Majorana zero-modes are 
localized near the ends of the wire (see also Fig. 12), while 



the rest of the states extend throughout the entire system. 
The corresponding LDOS is shown in Fig. [29] 



surements should suffice in establishing the existence of 
the zero-energy Majorana edge modes in semiconductor 
nanowires. 



A. Local density of states in superconducting 
nanowires 

In the previous section we have shown that disorder 
generates low-energy states and reduces the mini-gap. 
Nonetheless, a small mini-gap does not implicitly mean 
that the Majorana bound state cannot be resolved in a 
spectroscopic measurement. The key observation is that 
the undesirable low-energy states are generally localized 
near impurities and defects (see, foe example. Fig. 21). 



A local measurement could easily distinguish between a 
zero-energy state localized at the end of the wire and a 
low-energy state localized somewhere inside the wire. To 
clarify the question regarding the real space distribution 
of the low-energy spectral weight, we calculate the local 
density of states (LDOS) for several relevant regimes and 
compare the LDOS of a clean ideal system with that of 
a disordered realistic nanowire. 

We start with a clean nanowire with three occupied 
sub-bands (/i = 54. ) coupled non-uniformly to an 



s-wave superconductor. The coupling parameter is char- 
acterized by a non-homogeneity factor ^ = 0.8 and an 
intermediate coupling strength 7 = Aq, where Aq is the 
SC order parameter of the superconductor. In the ab- 
sence of a Zeeman field (F = 0), the nanowire is a triv- 
ial superconductor. Increasing the Zeeman field F above 
a critical value induces a transition from the trivial SC 
phase with = (no Majorana modes) to a topologi- 
cal SC phase with TV = 1 (one pair of Majorana modes). 
Further increasing F drives the system though a series of 
alternating phases with trivial {N even) and non-trivial 
(TV odd) topologies (see the phase diagram in Fig. [9|. Of 
major practical interest are the first two phases {N = 
and N = 1), as stronger Zeeman fields involve smaller 
gaps (see Fig. [To| or may destroy superconductivity 
altogether. Typical spectra from the first three phases 
(A^ = 0, 1, 2) are shown in Fig. 28 and the corresponding 
LDOS is shown in Fig. 
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The main conclusion suggested by the results shown 
in Fig. [29] is that clear-cut evidence for the existence of 
the Majorana zero modes can be obtained by driving the 
system from a trivial SC phase with A/" = to a topo- 
logical SC state with = 1 by tuning the Zeeman field. 
In the trivial SC phase there is a well-defined gap for 
all excitations, including states localized near the ends 
of the wire. By contrast, the topological SC phase is 
characterized by sharp zero-energy peaks localized near 
the ends of the wire and separated from all other exci- 
tations (including possible localized in-gap states) by a 
well-defined mini-gap. 

Is it possible to clearly distinguish the two phases with 
different topologies in the presence of disorder? The an- 
swer is provided by the results shown in Fig. 
a nanowire with charged impurities. In contrast 



30] for 
with 

the clean case, all the low-energy states are strongly lo- 
calized. Nonetheless, the signature features of the two 
phases (the gaps and the zero-energy peaks) are pre- 
served. At this point we emphasize two critical proper- 
ties: i) The features illustrated in Fig. 30 are generic, i.e., 
they do not depend on the type or the source of disorder. 
Similar LDOS can be generated using any other signifi- 
cant type of disorder discussed in the previous section, or 
combinations of different types of disorder, ii) Observing 
a zero-energy peak at a certain value of the Zeeman field 
does not by itself prove the realization of a topological 
SC phase. The trivial SC state with N = 2 may also 
have a zero-energy peak separated from all other excita- 
tions by a mini-gap. To clearly identify the A^ = 1 phase 
one must measure the LDOS as a function of the Zeeman 
field starting from F = 0, i.e., from the trivial SC phase 
with A^ = 0. Continuously increasing F will generate a 
transition from a phase characterized by a well-defined 
gap to a phase with strong zero-energy peaks localized 
near the ends of the wire. But what is the signature of 
the transition? 



Figure [31] shows the spectrum and the LDOS of a sys- 
tem with a Zeeman field F = 21.5£^q,. In a clean wire 
this corresponds to the transition between the N = and 
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FIG. 29: (Color online) LDOS for a clean nanowire in three 
different phases: trivial SC phase with A/" = (F = 15^a, 
top), topological SC phase with = 1 (F = 25Ea, middle), 
and trivial SC phase with AT = 2 (F = 45^^, bottom). The 
corresponding spectra are shown in Fig. |28| Notice the finite 
energy in-gap states localized near the ends of the wire (top 
and bottom) and the zero-energy Majorana modes (middle 
and bottom). The weight of the zero-energy modes in the 
N — 2 phase (bottom) is twice the weight of the Majorana 
modes in the topological SC phase with = 1 (middle). 
However, the clearest distinction can be made between the 
N — and the N — 1 phases. The LDOS is integrated over 
the transverse coordinates y and z. 




FIG. 30: (Color online) LDOS for a nanowire with charged 
impurities. The top picture corresponds to a trivial SC state 
with N — 1^ while the bottom picture is for a system with 
N — 1. The linear impurity density is riimp = 7 / /jm. Notice 
that all the low-energy states are strongly localized, but the 
clear-cut distinction between the two phases holds. 

type of critical behavior will characterize a finite range 
of Zeeman fields. Observing the transition between the 
topologically trivial and nontrivial phases, which is char- 
acterized by the closing of the gap and by a spectral 
weight distributed over a wide energy range, is the final 
ingredient necessary for unambiguously identifying the 
Majorana bound states using LDOS measurements. 



= 1 phases, which is marked by the vanishing of the 
quasi-particle gap, as shown in figures [To] and [121 The 



LDOS is characterized by a distribution of the spectral 
weight over the entire low-energy range of interest. This 
property holds at any position along the wire. Adding 
disorder induces localization, but does not change this 
key property. In fact, based on the analysis of the results 
shown in Fig. [241 we know that in disordered systems this 



B. Tunneling differential conductance 

An ideal type of measurement that exploits the proper- 
ties identified in the previous subsection consists of tun- 
neling into the ends of the wire and measuring the dif- 
ferential conductance [72H76] . To a first approximation, 
dl/dV is proportional to the local density of states at the 
end of the wire, so the general discussion presented above 
should apply. Here, we focus on certain specific aspects 
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of a tunneling experiment, e.g., the specific form of the 
tunneling matrix elements and the role of finite tempera- 
ture, that may limit the applicability of our conclusions. 
We find that values of the parameters consistent with an 
unambiguous identification of the Major ana bound state 
in the semiconductor nanowire are well within a realistic 
parameter regime. 



FIG. 31: (Color online) Energy spectrum and LDOS in the 
vicinity a topological phase transition. The top panel shows 
the spectra of a system with F = 21. 5 without disorder 
(blue diamonds) and in the presence of charged impurities 
(riimp = 0.7//im, orange squares). Note the absence of a 
gap. The corresponding LDOS are shown in the middle (clean 
system) and bottom (disordered wire) panels. The spectral 
weight is distributed over the entire energy range, including 
at positions near the ends of the wire. 



The tunneling current between a metallic tip and 
the nanowire can be evaluated within the Keldysh non- 
equilibrium formalism [77]. In terms of real space Green's 
functions we have 



I=^j douReTr { [l - G^{uj)T''{uj - eV)] ' [(1 - 2/^_ey)G^(c^)r^(a; - eV) (46) 

+2(^_ey - fu)G^r^{uj - eV) - (1 - 2U)G^{uj)T^{uj - eV)] [l - G^{uj)T^{uj - eV)] , 

I 

where V is the bias voltage applied between the tip and the nanowire and = l/(e^^ + 1) is the Fermi- 
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FIG. 32: (Color online) Differential conductance for tunneling 
into the end of a superconducting nanowire. The curves cor- 
respond to different values of the Zeeman field ranging from 
r = llEc, (bottom) to r = 36^a (top) in steps of Ec,. The 
curves were shifted vertically for clarity. The trivial SC phase 
(r < 21^0,) is characterized by a gap that vanishes in the criti- 
cal region (F ^ 21 E^)- The signature of the topological phase 
is the zero-energy peak resulting from tunneling into the Ma- 
jorana mode. The differential conductance was calculated at 
a temperature T ^ 50mK for a disordered wire with a linear 
density of charged impurities — 1 1 \ivn. 



Dirac distribution function corresponding to a tempera- 
ture /c^T = The retarded (advanced) Green's func- 
tion for the nanowire has the expression 



r ,uj 



En±ir] 



(47) 

where Un and Vn are the particle and hole components of 
the wave function corresponding to the energy En. The 
wave functions and the energies are obtained by diago- 
nalizing the effective BdG Hamiltonian for the nanowire, 
including the contibution from disorder, as described in 
the previous sections. The matrices F^^^^ contain in- 
formation about the tip and the tip-nanowire coupling. 
Specifically, we have 



F^(^)(r,r^c^) 



dx- 



(48) 



X ±ir] 

where u{x) represents the density of states of the metallic 



tip and jr,r depends on the tunneling matrix elements 
between the tip and the wire. We note that in Eq. (46) 



the trace is taken over the position vectors. We consider 
a tunneling model in which the amplitude of the tunnel- 
ing matrix elements vary exponentially with the distance 
from the metallic tip. Specifically, we have 7r,r7o^r^r', 
where 70 gives the overall strength of the tip-nanowire 
coupling and the position-dependent factor is 



"I ix-xtipy-\-iy-ytipy-\-iz-ztip)'^-xtipj 



(49) 



with {xtip^ytip^ Ztip) being the position vector for the tip 
and ^ a characteristic length scale associated with the ex- 
ponential decay of the tip-wire coupling. In the numer- 
ical calculations we take ^ = 0.4a and {xup^yup^ zup) = 
(— 3a, L2/2), i.e., the tip is is at a distance equal 
with three lattice spacings away from the end of the wire. 
With these choices, the differential conductance becomes 



dV 



■ ^ [fiEn - eV)\{un\0)\' + n-En - eV)\{vn\0)\' 

n 

(50) 

where the matrix elements {un\0) and {vn\0) involve sum- 
mations over the lattice sites of the nanowire system and 
provide the amplitudes for tunneling into specific states. 
Finite temperature effects are incorporated through the 
derivatives of the Fermi-Dirac function, /^ 

The dependence of the tunneling differential conduc- 
tance on the bias voltage for a superconducting nanowire 
with disorder is shown in Fig. 32 Different curves cor- 



respond to different values of the Zeeman field between 
F = llEa (bottom) and F = 36Ea (top) and are shifted 
vertically for clarity. The temperature used in the cal- 
culation is bOmK, a value that can be easily reached 
experimentally. Lower temperature values will generate 
sharper features, but the overall picture remains quali- 
tatively the same. Note that the closing of the gap in 
the critical region between the trivial SC phase and the 
topological SC phase can be clearly observed. In this re- 
gion dl/dV has features over the entire low-energy range, 
as discussed in the previous subsection. The Majorana 
bound state at F > 22Ea is clearly marked by a sharp 
peak at y = 0, separated by a gap from other finite 
energy features. We conclude that measuring the tun- 
neling differential conductance can provide a clear and 
unambiguous probe for Majorana bound states in semi- 
conductor nanowires. 



VI. CONCLUSIONS 

In conclusion, we have developed a comprehensive the- 
ory for the realization and the observation of the emer- 
gent non-Abelian Majorana mode in semiconductor (e.g., 
InAs, InSb) nanowires proximity coupled to an ordinary 
s-wave superconductor (e.g., Al, Nb) in the presence of a 
Zeeman splitting induced by an external magnetic field. 
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The importance of our work lies in the thorough investi- 
gation of the experimental parameter space, which is re- 
quired in order to predict the optimal parameter regime 
to search for the Majorana mode in nanowires. Since the 
number of possible physical parameters in the problem 
is large (e.g., electron density in the nanowire, geomet- 
ric size of the wire, chemical potential, the strength of 
spin-orbit coupling, the superconducting gap, the hop- 
ping matrix elements between the semiconductor and 
the superconductor, the strengths of various disorder in 
the semiconductor, or in the superconductor, or at the 
interface between them), theoretical guidance, as pro- 
vided in this work, is highly desirable for the success 
of the experimental search for the Majorana fermion in 
solid-state systems. Aside from the obvious conclusions 
(e.g., strong spin-orbit coupling and large Lande ^-factor 
in the semiconductor, large superconducting gap in the 
superconductor, and low disorder everywhere stabilize 
the topological phase), we have discovered several unex- 
pected results. In particular, we find somewhat surpris- 
ingly that the strict one-dimensional limit with purely 
one-subband occupancy for the nanowire, as originally 
envisioned by Kitaev [TT] and later used by many re- 
searchers [m [371 EH], is not only unnecessary, but is 
in fact detrimental to creating Majorana modes. In the 
presence of disorder, the optimal system should have a 
few ( 3-5) occupied subbands in the nanowire for the 
creation of the Majorana modes with maximal stability. 
This is, of course, great news from the experimental per- 
spective, because fabricating strictly ID semiconductor 
nanowires with pure one subband occupancy is a chal- 
lenging task. Another important result of our analysis 
is the relative immunity of the Majorana modes to the 
presence of disorder in the system. The most dangerous 
disorder mechanisms are due to charged impurity cen- 
ters in the semiconductor and to inhomogeneous hopping 
across the semiconductor-superconductor interface. Our 
work suggests that the optimal nanowires for observing 
the Majorana mode should not only have as little charge 
impurity disorder as possible in the semiconductor, but 
they should also have a thin insulating layer separat- 
ing the semiconductor and superconductor (so that the 
tunneling across the interface is not too strong), as well 
as some non-uniformity in the tunneling amplitude be- 
tween the semiconductor and the superconductor across 
the width (but not the length) of the nanowire. Our 
detailed numerical calculations establish that the zero- 
energy Majorana modes should clearly show up in exper- 
iments, even in the presence of considerable disorder in 
the nanowires. In addition, we establish that the disor- 
der in the superconductor has little effect on the Majo- 
rana mode in the nanowire. Another salient aspect of our 
work is the detailed calculation of the expected tunneling 
spectroscopy spectra for observing the Majorana mode 
in the nanowire using realistic physical parameters. Our 



results establish that the predicted topological quantum 
phase transition between the trivial phase with no Ma- 
jorana mode to the topological phase with a well-defined 
zero energy Majorana bound state should be clearly ob- 
servable as a striking zero-bias anomaly in the tunnel- 
ing current when the Zeeman splitting is tuned through 
the quantum critical point separating the two phases. 
More importantly, we calculate realistic tunneling spec- 
tra in the presence of uncontrolled spurious bound states 
in the system which are invariably present in real samples 
due to the localized random impurities in the semicon- 
ductor environment, clearly showing how to discern the 
topological features associated with the Majorana bound 
states from the background of contributions due to trivial 
bound states caused by impurities in the system. 

Our work emphasizes the tunneling measurements, 
which would directly establish the existence of a ro- 
bust zero energy mode in the system, providing the 
necessary condition for the existence of the Majorana 
fermion. What we have done here is to develop a detailed 
theory for the existence of a topological phase in the 
semiconductor-superconductor heterostructure, by tak- 
ing into account essentially all of the relevant physical 
effects. Once the presence of a robust zero energy mode 
is established, hence the necessary condition for the ex- 
istence of the Majorana is realized, one must move on 
to establish the sufficient condition, which would obvi- 
ously be a harder task. Several ideas for establishing 
definitively the existence of the Majorana mode (and its 
non-Abelian braiding statistics nature) have already been 
suggested in the literature, including experiments involv- 
ing the fractional Josephson effect [11, 14j|33l[37l[78], the 
quantized differential conductance [76^.79], and Majorana 
interferometry [8QH82] . 

Our work establishes the realistic likelihood of the ex- 
istence within laboratory conditions of the non-Abelian 
Majorana zero-energy mode in spin-orbit interacting 
semiconductor nanowires proximity-coupled to ordinary 
superconductors. We also establish that tunneling spec- 
troscopy is one of the easiest techniques to directly ob- 
serve the elusive Majorana in realistic solid-state sys- 
tems. Greater challenges, such as carrying out topologi- 
cal quantum computation, lie ahead once the laboratory 
existence of the Majorana mode is established experimen- 
tally. 
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